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Electromagnetic Field Theory Solution of the Infinite 
Tapered-Plane Transmission Line 


By Noacu Amiray, ABRAHIM Lavi, and FREDERICK YOUNG, Pittsburgh, 
Pennsylvania, U.S.A.1) 


Introduction 


Aside from a few isolated cases an exact solution to the nonuniform tapered 
transmission line problem is very difficult and consequently certain simplifying 
approximations become mandatory in arriving at a reasonable accurate solu- 
tion. The methods employed in the literature are based mainly upon the tele- 
graphy equations. The series and shunt impedances are obtained either from 
the step-line approximation where the line is broken into segments of parallel 
planes, or are assumed to vary according to an arbitrary function. The cor- 
relation between the mathematical solution and the actual physical structure is 

lacking in many cases. In this work a new approach is presented which bridges 
the missing link between the circuit theory analysis and the line geometry. The 
mathematical technique used is based upon the exact electromagnetic field 
theory solution of the oblique-planes line in conjunction with LAGRANGE’s 
method of variation of parameters. 


Wave Propagation between Two Infinite Oblique Conducting Planes 


Prior to the analysis of an arbitrary tapered-plane transmission line, the 
solution of the two infinite oblique-planes line of figure 1 is briefly stated. 
Assuming an isotropic and nondissipative medium, ¢ and uw are both constant 
scalars, and assuming the planes to have a very high conductivity, the solution 
of the TM modes of propagation is given by 


H, = [C,,, HD (kr) + Cy, HO (k1r)] (A,, cosm6 + B,, sinm 6) , (1) 


E=—"—_[C,,, Her) + C,,, Hk r)] (B,, cosm60 —A,,sinm6), (2) 


z. 1MeEerv 


14 CAM (kr) + Cy, H(kr)] (A, cosmO + B,, sinm6), (3) 


9 iwe ad 
where H!)(k vr) and H'?)(k 7) are Hankel functions of the first and second kind 


1) Carnegie Institute of Technology, Department of Electrical Engineering. This work was 
submitted by Noacn Amiray in partial fulfillment of the Degree Doctor of Philosophy. 
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2m? 


of mth order and Cym; Cami Ami Bm and m are constants of integration and | 
separation. k? = w? we = (27/A)?, A being the free space wave length. | 


H 
1 


Figure 1 
Two infinite oblique planes. 


For a TEM mode m = 0 and hence 
H,=C, Hk + CA (en, (4) 
E,=—iV/# [C, HMR + C, HRD). (5) 


It can easily be shown that higher order modes will be eliminated if 7 0) < 4/2. 


Wave Propagation in an Infinite Tapered-plane Transmission Line 


Consider the case of two conducting sheets extending from z= —oo to 
z= +00 as shown in figure 2. One plane lies along the x-axis (y = 0) while 


yh 


0 Q@ xX X+AXx x 
Figure 2 


Infinite tapered plane transmission line. 


the other is given by y, = f(x) where f(x) is an arbitrary function of x. Draw 
tangents to the tapered plane y, at the points M and N of figure 2. These 
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tangents intersect the x-axis at O and Q, while intersecting one another at the 
. point F. Designate the length OM as v,, the length ON a8 7,44, the angle 
between OM and the x-axis as 6, and the angle between QN and the x-axis as 
fe... ben 


_ WN dy 
tan 6, = Fr and tan 0,4 0% = Ge les ay (6) 
Divide the tapered line into small sections VF, FG and so on, as shown in 
figure 3. Each section is a portion of an oblique plane making an angle of 0, with 
the x-axis. The electromagnetic field quantities in each section are given by 


(4) and (5). 


0 Q X -X+AX x 
Figure 3 


Two small sections of an infinite tapered-plane transmission line. 


It should be emphasized here that 7, and 6, are functions of x for a certain 
y, = f(x) and thus the coefficients C, and C, of (4) and (5) vary from one 
section to another depending on x while being constants over each section. 
The field quantities everywhere between the planes are determined by requiring 
the magnitudes of the corresponding electromagnetic field components at the 
end of a section to equal those at the beginning of the following section, while 
taking the limit as Ax > 0. 


Consider the section VF of figure 3. With center O and radius OF draw the 
arc FT, where T is its point of intersection with the x-axis; with center Q and 


radius OF draw the arc OL, where L is its point of intersection with the x-axis. 
Since the mode of propagation considered is the TEM, the arcs FT and FL 


represent the wave fronts for the sections bounded by VF and FG respectively. 
The corresponding electromagnetic field quantities are matched along these 
arcs as stated above. | 

In a typical section VF, 


H,, = C,(x) Hy (Rr) + Cy(x) A(R 7) , (7) 
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In the adjacent section FG, 


H,, = C,(% + Ax) HD) (kr) + C(x + Ax) He (kr), (9). 


z 


Ey. = —i Z[C,(x + Ax) HM (kr) + C(x + Ax) HY (R7)], (10) | 


where OF 7S Ce 
For values of ky > 1 or g, the Hankel functions can be replaced by their 
asymptotic form 


H (kr) = (<5 -)" gilkr—2a+1)a/4] (11) 
2 1/2 —if[kr— +1)z, 
Hk r) = (<5 -) e-iler—(2q+1) 2/4] | (12) 


For g = 1 and kv = x the error introduced is about 3%. 
Matching the corresponding field quantities along the arcs FT and FL 
while Ax > 0 implies 


jim Fal, -oF = jim Hal, oF i Fe 
Jim Foil, -o8 = jim Fool, -oF ; (14) 
From figures 2 and 3 
eae a Ax 
OF = 1, + cos6, ’ ee 
+ _ (= wae 
QF Ve 4Ax cos aes , (16) 


where 0 Xa <1. 


Upon the substitution of the asymptotic form of the Hankel functions (13) 
and (14) become 


(l—a) Ax a 


= lim hd eae [lee Fay eps ese carsosyer 
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2 aAx 38a 
: Ki ley oi [Met esee) “| 
zt kR Y ae a Ax ilk adAx 32 
2 cos@ a (re cos 6 | 
x + C(x) e x 


__(1=a) Ax )-22] 


aN AEM Ye a aN 


fy 


4 


(1—a) Ax _ 38a 


+ C,(*% + Ax) a [Rrra el at : 


The mathematical solution is carried out by expanding C,(x + 4x); 
_ C,(% + Ax); 7,44, and 6,,4, into Taylor series around the point x. Taking the 
_ limit as Ax +> 0, (17) and (18) produce two differential equations for the coeffi- 
cients C,(x) and C,(x) given by 


28) C1) (Ft - Az) (F-- +), (19 
Seen let) 


Integration of (19) and (20) yields 


dx : dx 
7, cos6, tag 


log C,(x) = log A + — + logr, —tkr,— 2 | 


logC,(x) = log Bar * logr, +tkr,— al ae —ikf okay (22) 


2 v,, cosé, 


where logA and logB are constants of integration. 
Referring to figures 2 and 3, the following relations can be derived. 


aaa (BY (BYE 


1 
cos?6, 


—1+ tan?6, =14+(42)", (24) 


The use of (23) and (24) in the evaluation of the last two integrals in (21) 
and (22) gives 


7 (2) le Vy ay =f t=} 
Y esi =/[+(@ We /(ai) b+ Ge] FE peas 


dx (25) 
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[ hp- f+ (GNP anne 9) 


Xo 


where s, , is the length of the tapered line arc between an arbitrary fixed 
point [x93 ¥1(%o)] and the point [x; y,(x)]. ' 

In this case, the choice is %) = 0. Thus, (21) and (22) assume the form which 
is given by 


Cy(x) = A ae oe *(8s-%s) , (27) 
V1 

C(x) = B eaig en tBu") (28) 
V1 


The final equations of the electromagnetic field quantities are 


nom (che) lalte)totlr tow «g(a ene tenal 


2 RY, Vi Vi 


= ee [A, elk sx BS A, gti , 
37 32 


E,= iz 2 if ( "x i lA flare! PB ed bea | | 


= 2° (— Ay et k Sx =a A, gf FBSy) . | 


(30) 


where A, = A (2/7 k)'* e—*#/4 and A, = B (2/a k)'? e*™* . 

To distinguish this solution from the exact field solution derived directly 
from MAXWELL’s equations, and from the solution obtained from the telegraphy 
equations, this solution will be henceforth designated as the ‘Approximate field 
solution’ of the tapered line. 

Since energy is conserved and the transmission line is assumed to be lossless, 
the instantaneous power flow per unit of z width through any transverse cross- 


section extending between the two planes should be constant. The RMS power 
flow is 
0 


Piss : fReLE, H*)\ 7 aG = const. (31) 
0 


Upon the substitution of (29) and (30) into (31) 


| (32) 
=F W(Aelt — 4.) 780 = 5 (Asi — 149) | 


ial 


Vol. XIT, 1961 Electromagnetic Field Theory Solution 95 


It appears that the power P flowing through any transverse cross-sectional 
“surface i is a function of x. The dependence of P on x expressed by y, = 6,/sin 6, 
violates the law of conservation of energy. 

However, y, is very close to unity for 0 < 0, <z/6 (at 0, = 2/6; is 1-04). 

Therefore, if 0, is kept below z/6 for any point along the ree no serious error 
i is introduced. 
_ This limitation on the tapering angle is not surprising if it is realized that 
in the procedure of matching the corresponding field quantities between two 
adjacent segments (see (17) and (18), the wave fronts were assumed to be 
_ cylindrical surfaces. Obviously this assumption breaks down at large angles 
_which are not common in the art. Hence, for all practical purposes the law of 
_ conservation of energy is not violated. 


aot yee eae \ yO AN rece Cee 


Transformations by the Nonuniform Line 


| Any transmission line that may be considered as a four terminal network 

has the property of impedance transformation. Thus, designating the ratio 

of input voltage to input current as the input impedance, Z;, and the ratio of 

output voltage to output current as the output impedance (terminating imped- 

ance), Z,, the impedance transformation ratio is given by Z,/Z;. 

. An the general transmission line of figure 4, the MES between points P, 
and P, is 


: ; 70, 
V, -/e, yd = Z (A, e~** — A, efhss) —*_* | (33) 
é Vv 


0 xf Xs; x 


Figure 4 


General transmission line. 


The current flowing in a strip of width d is 


Aer ae Ae") (34) 
V1 


For a finite length line having the total arc length 


GS -/+(§ (22)" [Pax 
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as in figure 4, A, and A, are expressed in terms of the voltage and the current . 
at the end of the line (V; and I, respectively). . 
Thus, the voltage and current at any point x along the line are 


v= g (22 - (4 Is me, x) gik(S—sy) _ (4 - Tr) eins 1, Oe 5 : 
(35) 


te F(R) (+ tig on + Ub ake) Ae). 


In order to suppress reflection in the line, the following relation must hold 


Z,- p= 2. bs (37) 


In other words, if the line is terminated with a load impedance Zs, no 
reflections occur. In this case the input voltage and current (x = 0) are 


Wav. S a gies E94 Ts (2) giks (38) 
Ve 05 me 7 eNaS : 
Te OVGT aks 2 es Vos ikS 
is iG Zr05 a, telnbes nee 


and the input impedance is 


| EO A PE 


: 40 
I; Is ys S ys be 
The impedance transformation ratio for this case becomes 
LS WS 
z=, 41 
Zi Sh (41) 
While the voltage and current transformation ratios are 
Ve. Sag \ He” as 
vom ine \ nares ie 
Tse a Ys \ Paes 
To ee 43) 


It is interesting to note that the transformation ratios are explicitly ex- 
pressed in terms of the function y,, and the impedance transformation ratio 
depends solely on the ratio of the output to input vertical separation between 
the lines. The actual shape of the line has no bearing whatsoever on the imped- 
ance transformation ratio; that is, given two lines described by y, = f,(x) and 
¥1 = fe(x) such that /,(0) = f,(0) and f,(L) = f,(L) as in figure 5, the impedance 


transformation ratio for both lines is identical as long as both lines satisfy the 
conditions imposed previously. 


“ANA Ree 


ek Oe laa Lah Rie 
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The previous conclusion has to be modified, however, when the voltage (or 
_ current) transformation ratio is evaluated. The absolute value of the ratio is 


; again independent of the function, y,, but there is definitely a phase ere 
upon the length of the arc S. Therefore, for the two lines f,(%) and f,(x) o 
efigure 5, if 


i 


fs (Ge yf aap GeyP ews 
0 0 


_ then the two lines do not behave identically phase-wise although /,(L) = /,(L). 


y 
Yah) 
L Xx 
Figure 5 


Two transmission lines having the same impedance transformation ratio. 


Design Considerations 


The limitation 0,,., <2/6 in conjunction with 7/A > 1/2 and 0,,,,7/A < 1/2 
define the region within which the solution is valid. These three conditions are 
represented in figure 6. The shaded area gives the permissible combinations for 
rv/A and 0,,,,- It is interesting to note that in order to increase the bandwidth 
while retaining the same impedance transformation 6,,,, must be decreased. 


0 10 20 30 40 
Gnax (degrees) 


Figure 6 
Impedance transformation region. 
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Comparison with other Existing Solutions 


The solutions of other tapered lines obtained from the telegraphy equations 
check in certain respects with those obtained from the approximate field” 
solution. It has been found that the impedance transformation in Bessel and 
exponential lines depends upon the ratio of the vertical separations between 
the electrodes at the two ends of the line. However, a discrepancy in phase is 
observed. The phase variation in the approximate field solution depends upon 
the length of the arc while in the telegraphy equations solution the phase 
variation depends upon the linear distance x. This error in phase is not sur- 
prising in view of the fact that for the simple case of the two oblique-planes 
line the phase variation obtained from the telegraphy equations solution is in 
disagreement with the solution obtained directly from MAXWELL’s equations. 
On the other hand in the approximate field solution, the phase variation is 
expected to be more accurate because the shape of the wave fronts is accounted 
for. Moreover, the boundary conditions (i. e. the electric field is perpendicular 
to the conducting surfaces) are satisfied in contrast with the situation in the 
step-line approximation. 

The agreements between the telegraphy equations solution and the approxi- 
mate field solution are encouraging but certainly not profound. While the 
approximate field solution yields quick and accurate results for any arbitrary 
tapering function y,, the telegraphy equations in conjunction with the step 
line approximation lack this generality. 


Conclusion 


A new method for the study of the steady state performance of infinite 
tapered-plane transmission lines has been presented. This method is advanta- 
geous because it offers a quick and accurate evaluation of the transformation 
properties of the line. The impedance transformation ratio as well as the 
absolute values of the voltage and current transformation ratios of a matched 
line depends upon the ratio of the initial and final separations regardless of the 
tapering function of the line. In contrast to previous works, this paper accounts 
for the phase variation which is proportional to the tapered-line arc. 

In addition, design criteria are presented, which limit the tapering of the 
line and determine the bandwidth. 

The results derived in this work hold true for tapered strip lines as well as 
coaxial cables where the radial separation between the electrodes is small 
compared to the radius of the electrodes. 
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Zusammenfassung 


Fiir den eingeschwungenen Zustand der tapered-plane transmission line wird eine 
neue Lésungsmethode beschrieben und gezeigt, dass die Impedanztransformation 
direkt vom Verhaltnis der Abstande zwischen Anfang und Ende der leitenden 
Ebenen bestimmt wird. Die Phasendifferenz der Spannungen zwischen den Lei- 
tungsenden ist dem Offnungswinkel der beiden Ebenen proportional. Kriterien fiir 
den Entwurf solcher Leitungen lassen sich angeben. : 


(Received: August 30, 1960.) 
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Duct Flow in Magnetohydrodynamics 


By Curen C. Canc and THomas S. LunpcREN, Minneapolis, Minnesota, USA?) | 


Introduction 


This paper is an extension of the work of HARTMANN [2] *) and SHERCLIFF 
(3, 4] on the steady flow of conducting fluids through ducts under transverse 
magnetic fields — the simplest class of magnetohydrodynamic problems. We 
are concerned here mainly with the boundary value problems associated with 
flow in ducts with conducting walls. 


Equations and Boundary Conditions 


The following set of vector equations appears to give an adequate description 
of the steady state interaction between electromagnetic and hydrodynamic 
forces: 


curl Bia ped (1) 
div B=0, (2) 
curl H=0, (3) 
div E=o,/é&, (4) 
jo= OBS x Bee (5) 
oV-VV=-Vp+orVV+j7xB+o,E (6) 
div, V = 0% (7) 


It is assumed that the magnetic and dielectric properties of the medium are the 
same as in a vacuum; fly and €) are the magnetic permeability and dielectric 
constant in vacuum. In (5), Oum’s law, the electrical conductivity o is assumed 
constant for a homogeneous medium. Equations (6) and (7) are the momentum 
and continuity equations which describe the steady motion of an incompressible 
fluid. If the system to be analyzed is composed partly of fluid and partly of 
solid or vacuum, the last two equations, (6) and (7), only have to be satisfied 
in the fluid, while the first five equations must be satisfied throughout all space. 
The description of the system is completed by specifying zero velocity at rigid 
boundaries and imposing continuity of tangential components of E and of 


1) Department of Aeronautical Engineering, University of Minnesota. 
*) Numbers in brackets refer to References, page 114. 
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-hormal and tangential components of B at interfaces, plus boundary conditions 
4 on B and E at infinity. It should be mentioned that continuity of tangential 
“components of B implies no surface currents. This is the case, since in non- 
“magnetic materials surface currents occur only in the presence of unsteady 
Magnetic fields. 

_ In equations (5) and (6) the free charge o, will be neglected. It will be 

retained in (4), which serves only to determine the free charge once E is known. 
This can be considered as the first step in an iterative process, the second step 

being to substitute the value of g, calculated in the first step into (5) and (6). 


0-0, 


Region 30=0 


Teg ad et dB 


Figure 1 
Cross-section of duct. 


Figure 1 is a sketch of the system under study. An electrically conducting 
incompressible fluid (region 1) flows in the z-direction through a straight duct 
whose walls are of constant thickness (region 2). The electrical conductivity 
of the wall is o,. Outside of the duct (region 3) the conductivity is zero, and at 
infinity a uniform magnetic field By acts in the y-direction. 

If it is assumed that the velocity has only a z-component, that all physical 
quantities (except pressure) are independent of z, and that there is no net flow 
of current in the z-direction, then it can be shown that B, = 0, B, = Bo, and 
pois =.0, 

Using the results of the previous paragraph, the vector equations from (1) 
to (7) can be reduced to two second order linear partial differential equations, 
namely, the z-component of the momentum equation 

0 B, OB, 
= 7 (8) 


e 


and the z-component of the curl of OHM’s a 


OV, 
V?B, +o Mo Be pe + (9) 
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In the first of these 0f/0z is a constant, since from (6) Vp is independent of z. 
In the fluid, region 1, both (8) and (9) are valid, while outside of the fluid, 
regions 2 and 3, only (9) is true with V, = 0. In addition, since region 3 is non- 
conducting, the current must be zero there. This fact and (1) imply that B, is 
constant. Therefore, since B, goes to zero at infinity, B, must be identically 
zero in region 3. 

Across the boundaries of the several regions, B, and the tangential compo- 
nent of electric field must be continuous. Using OHm’s law, the condition on 
the electric field can be written 

[|= 


where 7, is the current in the direction of the tangent to the interface. (A 
square bracket around a quantity means the discontinuity in this quantity.) 
This is true since the velocity vanishes at the boundary, causing the V x B 
term in OHM’s law to vanish. Observe now that the x- and y-components of 
(1) are 

‘ OB, ; OB, 
Mo Ix = “py > Mo ly = — Ox 


so B, can be considered as a ‘stream function’ for “oj — the current flows 
along the lines B, = constant. This shows that 


OB, 
Ho Ji a on? 


where n is the outward normal. Using this result, continuity of E, becomes 


Before summing up the results of the last few paragraphs, let new dimen- 
sionless variables be introduced by 


E eT 
| aa 2 ’ 
es Vz 
i (a*/» @) (dp/dz) ’ (10) 
B=— B, 


(a?/v @) (0p/0z) wy (veo) ’ 
ne o \1/2 
Melee 


oe 
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where a is a characteristic dimension of the duct and M is the Hartmann 
number. With this new notation and using subscripts 1 and 2 to denote the 
-tegion, the problem is to solve 


, 


; OB 
: V-V,+M os =—1, (11) 
. av, 
7 V?B,+M a = (0, (12) 
‘in region 1, and 

V?B, =0 (13) 


in region 2. The boundary conditions are: B, = 0 on C,, the boundary between 
2 and 3; and V;=0, B, = B,, o, 0B,/0n =o, OB,/0n on C,, the boundary 
between 1 and 2. 

This problem is difficult since it involves two domains and two sets of 
boundary conditions. In special cases it can be simplified. If the duct wall is 
a perfect insulator (0, = 0), then B, = 0 so that it is only necessary to solve 
(11) and (12) with V,=0 and B,=0 on C,. If the duct wall is a perfect 
conductor (0, = oo) the boundary conditions become: V, = 0 and 0B,/dn = 0 
on C,. There is another limiting case discovered by SHERCLIFF [4], for which 
the problem reduces to solving (11) and (12) in region 1 with boundary condi- 
tions given on C,. Suppose the thickness of the duct wall (h/a in the new nota- 
tion) is much smaller than unity. To a good approximation the harmonic 
function B, is locally linear, that is, it varies linearly across the duct wall. This 
can be seen clearly by considering the membrane analogy for solutions of 
LAPLACE’S equation. Then 0B,/0n = — a B,/h on C,. The boundary conditions 
become V, = 0, 0B,/0n + B,/p = 0, where py = o, h/o, a. It should be noted 
that the cases o, = 0 and og = oo are included in the last boundary condition 
with g = 0 and » = oo respectively. 


Parallel Sided Duct 


There is one situation in which the approximate boundary condition be- 
comes exact. This is the case of flow in a rectangular duct when the walls 
parallel to the applied magnetic field are at infinity. In this problem the har- 
monic function B, is independent of &, and hence must be a linear function of 7. 

With the dimensionless variables defined in the previous section, the pro- 
blem is to solve 


Vv dB aB ee 
pm cag May 


with boundary conditions V = 0, dB/dyn + B/y = 0 when y = + 1 respectively. 
Here, half the duct height has been taken as the characteristic length, a. The 
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solution of the above equations is easily found to be 


ape! gt+1 cosh M | 
ie iar M p + tanhM (t= cosh M ), | 


gt+1 sinh M 


Ei pul if Mo+tanhM coshM * 


(15) 


HARTMANN [2] gave this solution for the case g = 0. For given q, the velocity, 
which is parabolic for M = 0, becomes flatter in the center as M increases. 
When M is very large it tends to be uniform except in a boundary layer of 
thickness the order of 1/M. Asymptotically, 


Ue ee 
Ve aes (16) 
This shows that V ~ 1/M for non-conducting walls while V ~ 1/M? for per- 
fectly conducting walls: increasing the wall conductivity decreases the velocity. 

Another quantity of interest is M dB/dn, the ratio of Lorentz force to the 
magnitude of the pressure gradient. This is also proportional to the current 


density. By an easy calculation, 


dB +1 cosh 
ieee Pp tee 
dn fa Mg+tanhM cosh M 


This is plotted in Figure 2 for M=3, p=1, h= 0-54, to give a typical 
example. The value of this quantity in the wall, 


( M — tanh M 
a? ae SANT, 


15 
7] 


Figure 2 
The current distribution across half of a parallel sided duct for M = 3, p=1h=05a, 


RWS Oe eee 
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is also plotted. From this figure it is seen that near the center of the duct the 
current flows to the left. This gives a Lorentz force which tends to help the 
viscous forces balance the pressure gradient. Near the wall the current flows to 
the right giving a Lorentz force which opposes the viscous forces, so that the 
viscous forces have to be larger in order that the pressure gradient be balanced. 
For large M@ 


dB 


except in a boundary layer of thickness the order of 1/M, and in the wall. This 
shows that for large M, the Lorentz force completely balances the pressure 
gradient except in a thin boundary layer along the wall. 

Stated differently, the current distribution tends to be uniform to the left 
and of a magnitude such as to make the Lorentz force balance the pressure 
gradient. Since the total current flow must be zero, part of the return current 
flows in the boundary layer and part in the wall. When the wall conductivity 
is larger a greater proportion of the return current flows in the wall, taking the 
path of least resistance. This shows that for fixed M (large), less current flows 
in the boundary layer when the wall conductivity is high, indicating that the 
Lorentz force opposing the viscous forces is less. Therefore the viscous forces 
must be smaller in order that the pressure gradient be balanced. Now if the 
viscous forces near the wall are smaller, the second derivatives of the velocity 
will be smaller; hence the velocity of the core will be smaller. This shows why 
the velocity becomes smaller when the wall conductivity increases. 

Consider the dimensionless mass flow 


és. _ Mass flow per unit width 
qe [Va (a4]») (— 09/02) 


g+1 M-—tauhM 
M* Mo+tanhMmM ° 


The reciprocal of Q is essentially the pressure gradient required to maintain a 
given mass flow. In Figure 3 


Q(M = 0) 1 M* Mo+tanhM@ 


Q ~ 3 g+1 M-—tanhM (17) 


is plotted versus M for various values of y. This can be interpreted as the ratio 
of the pressure gradient to the pressure gradient required to maintain a non- 
magnetic flow with the same mass flow. Notice that for a given mass flow a 
much larger pressure gradient is required to maintain flow through a perfectly 
conducting duct than through a non-conducting duct. 


- 
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Figure 3 


The pressure gradient versus M for various values of p in a parallel sided duct. 


Perfectly Conducting Rectangular Duct 


In this section, the flow through a finite rectangular duct with perfectly 


conducting walls is considered. Regrettably, we have been unable to solve the. 


problem for arbitrary wall conductivity. 
The problem to be solved is 


OB OV 
2 et Mice 2 Eh a 
YeV+ M - iS a 0 


in a rectangle of width 2/ and height 2, with V = 0 on all the walls, 0B/d& = 0 
on the vertical walls, and 0B/0n = 0 on the horizontal walls. The expansions 


V= Doo, cosB, Hn, B=) b,sinBya, 1 =)" a, cosp,.4s 
7=0 j—0 =a 


where $,; = (7 + 1/2) 2, satisfies the boundary conditions on the horizontal 
walls and allows the differential equations to be written as ordinary differential 
equations for v; and b, — 

adv; ab; 


qe 7 ut MB; ava, Ge ger — BF by — M Bp, = 0. 


Solutions of these which satisfy the boundary conditions on the vertical walls are 


seek a; (1 7,3 Sinhy,;/ coshr,; sinhé + 7,,; sinhr,,; 1 coshr,; & 
J eke 2; COSh7,; 1 sinh7,; 1 + 7,; coshr,; 1 sinh7,; 1 


(18) 
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and 
i att a; ee topes sinhy,;/coshr,;& — v,; sinhr,; 1 coshr,, é 19 
J B? ae AU WN er, v,; coshr,; 7 sinhr,;/ 4- 7,; coshr,, 1 sinh, ; (19) 


j= (Bp +i MB), 75; = (6-1 MB) 


are complex. 7,; and 7,; may be separated into their real and imaginary parts, 
namely, 
rigs Ajctat 7; 
and 
19, =H —1Y;, 


where 
a; = (2)" [B; + (6? + MFR y= eat [— B; + (Bj) + M2epe, 


After some algebra to express v; and b; in terms of the real quantities «; and 
y; the final result becomes 


— yr 2 (= cosbyn fy _ 4s E,l8) — 75 Fil) 

ees B; B? + M? {1 a; sinh 2 a ei ane mit (20) 
_ wo 2(-1)) sinBjy {M a, F,(&) + ; E,(é) 

i Py By BR + M® tz, 3 a Sinh?2 aj 1 — 7, sin 2 If? 21) 


where 


E,(é) = “ [cosy, (J — &) sinha, (J + &) + cosy; (1 + €) sinha, (/ — €)], 
F,(é) = = [siny, (0 — &) cosha, (I + £) + siny, (1 + é) cosha, ( — €)]. 


The calculation of mass flow per unit pressure gradient is more easily 
accomplished by integrating the complex form of the velocity. The result is 
es il B; cosh2a;1—cos2a;/ | 


0= 81d (7 + M?) ( 1 (BP + M?)*? a; sinh2 a; 7 — yj sin2 yj J” 


(22) 


)(M = 0)/Q(M) is plotted versus M in Figure 4 for a square duct. The corre- 
sponding result for non-conducting walls, SHERCLIFF [3], is included for com- 
)arison. 
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Figure 4 


Pressure gradient versus M for mg = 0 and m = on, square duct. 


Flow in Arbitrary Symmetrical Duct for Large M 


‘It has been seen for rectangular ducts that the velocity distribution for 
large M consists of a uniform core with a boundary layer near the walls in 
which the velocity changes rapidly. The purpose of the present section is to 
investigate how the cross-section form and wall conductivity affect the velocity 
distribution in the core. This question has been answered for non-conducting 
walls by SHERCLIFF [3] and for circular ducts with walls of small conductivity 
by SHERCLIFF [4]. 

Let the duct be as indicated in Figure 1, except symmetrical about the 
x-axis. Let the upper surface be described by 7 = Y(&) the lower surface by 
ny = —Y(é). The problem is to solve 


OB 
POV Mime (23) 
OV 


with V=0 and 0B/0n+ B/p=0 at the wall. LetZ,=M(V+B),Z,=M(V—B). 
Adding and subtracting (23) and (24) gives 


: (WA 

ViZ pM he M, (25) 
OZ. 

Vad, — MM. (26) 


Singular perturbation theory, Levinson [1], says that at interior points and 
at points on the upper surface 


4, = Ly; + O(a) , 


ty 
ee 


wa 
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where Z,; is the solution of 
; OZ4; 
én 
which takes the given boundary values on the upper surface. There is a bound- 
ary layer near the lower surface and ends. Specifically, this is true when Z, 
is given on the boundary. In the present problem the value of Z, on the bound- 
ary is unknown but the above result should still be true, since if the exact 
solution of (23) and (24) were known the values which Z, takes on the boundary 
could be calculated and the problem formulated as a first boundary value 
problem. The analysis of Z, is similar, with Z, tending to the solution 0Z,, /On=1 
which assumes the boundary values on the lower surface. This reasoning gives 
asymptotic solutions 
Li; = Y(&) = Z,(E, Y(&)) , (27) 


Zo,;=H+ Y(E)+ Z(81, —Y(é)) , (28) 
where Z,(&, Y(é)) and Z,(&, —Y(é)) are the values which the exact solutions 


of (25) and (26) assume on the boundaries. Since Z,= M (V + B), Z,=M (V—B) 
and V is zero on the boundaries, it must be the case that 


Z,(&, Y(é)) = M BE, Y(), : (29) 

Z,(é, -Y() = — M Bit, -Y(®). 
But by the symmetry of the boundaries, B is an odd function of 7, therefore 
Z(é, -Y()) = M BE, Y(®). (30) 


Also, the velocity and induced field in the interior tend to 


»+ Zo; Y(é 
Pee Reve) +e (31) 
sae Y 
Bie = a va (32) 


In order to find B(&, Y(é)) it is necessary to have a relation between 
quantities across the boundary layer. By integrating V?B + M 0V/0n = 0 over 
a small cylinder which extends through the boundary layer, and using GREEN’s 
third identity, one finds 

ee a =) cos\n,.) <0, (33) 


On On 


| 
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where 1 is the outward normal. This result is equivalent to imposing continuity 
of tangential components of the electric field across the boundary layer. Now 
using V = 0 and 0B/0n + B/p = 0 on the boundaries, it is found that 


~ ee ne — M, cos (n, n) (34) 
on the boundary. Observe that 
PBs. OB cos(n, 8) + “i cos (x, 4) = — M-tcos(n, 7) (35) 
whence 
= = ([M-1— (M B(E, Y(€)) + ¥(&))] cos(n, n) (36) 


on the boundary. 

In particular this is true on the upper surface where B = B(é, Y(&)) and 
cos (”,y) = (1+ Y"?)-1/2. Substituting these into (36) and solving for B(é, Y(é)) 
gives 


M-1— Y 
BEE, Y(é)) aotae pee pt (i+ Y’)12 ° (37) 
With (37) and (31), V; is solved — 
_1 f ¥toatyy-in 
es eT ee ey: (38) 


This checks with (16) for the case Y = 1. Note, when the walls are non-con- 
ducting (y = 0) that V; = Y(&)/M as shown by SHERCLIFF [3]. In this case the 
velocity distribution has the same shape as the cross-section of the duct. In 
fact this result is also true for non-symmetrical ducts with non-conducting 
walls, that is, if the upper surface is 7 = Y,(€) and the lower surface is 
n = —Y,(6) then V; = (Y,() + Y.(&))/2 M. On the other hand for perfectly 
conducting walls (y = 00), V; = 1/M? which shows the velocity to be uniform 
in the core. Also, in this case, the dimensional velocity V, is independent of 
viscosity — V, = — (0p/0z)/(o B?). 
For a circular cross-section, Y = (1 — &2)1/2, a small calculation shows 


i+e  @-se 


Les M 1+4+Mo(i-®&i?° 


This differs from SHERCLIFF [4] in the occurence of the factor 1 + y; SHERCLIFF 
restricted y to be small. The volume flow through the circular duct is given by 
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1-2 
i+ Moa —pe & 


cos? 6 db (39) 
M 1+ M cos 6 
6 


“3,2 [4A 1 re m2 1 2 
M (M ¢)? (M ¢)* (M9)? 14+ Mo 


1—M 9\'/? 1— M \}1/2 
Sa ieee ee Ue ee eee 
ee sce Gear ‘ 


4 Note that the expression for Q, (39), does not have a singularity at M y = 0, 
the inverse powers are absorbed by the last term on the right. In fact when 
_ M gis small 


1+ Z 3 8 
Cee ig ee tas (MoS. 


In Figure 5 (Q/z) [4 M/(1 + ¢)] is plotted versus M ¢. Q/z is the average (dimen- 


sionléss) velocity. In Figure 6 V;/V,, is plotted versus € with M ~ as parameter. 
The latter figure shows the effect of wall conductivity on the shape of the 


07 
06 
05 


04 


No —— 
Figure 5 
4 M Vq,/(1+q@) versus M for circular duct, large M. 
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| 
velocity profile for the case of a circular duct. It should be noted in this case — 
that the shape depends only on the product M 9. 


—— 
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Figure 6 


Velocity distribution in a circular duct with M @ as parameter. 


The Effect of the Free Charge on the Electric Field 


In the introduction we stated that the free charge would be neglected in 
Ouw’s law and in the momentum equation, but would be retained in (4). We 
have seen that B, and V, can be determined without using (4), therefore E can 
be determined from Oxnm’s law leaving (4) for the determination of @,. This is 
done simply by taking the divergence of (5), with the result 


0, = — &) div (VV x B). (40) 

We shall now check to see that o, was negligible in (5) and (6) for the class 

of problems which we have considered. We note that in this case (40) becomes 
OV, 

Oc = Ey By “Aye (41) 


For large M it is apparent that @, will be largest in the boundary layer. To be 
more specific it will be largest near that part of the boundary which is parallel 
to the applied magnetic field. SHERCLIFF [3] has shown that the boundary layer 


has thickness of the order of alVM on walls parallel to the applied magnetic 
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- field. This allows us to estimate 0V,/0x by V, V,4,VM, ja at this point (V, ,, is 
the average or core velocity) which makes 


Z av a 


Pi iti ph ha eS) RTS bee 


0,— Of) By Vea I |. a ee 


“We note that in equations (5) and (6) the terms with @, will be small compared 
_tooV x B and 7 x B respectively provided the inequality 0,/o B, <1 is 
Z Ba tistied. With (42) this inequality can be put in the form 


; V, AR 
é - 7280 < c ? (43) 
4 e VM 


_ where c is the speed of light, A = wy ov and R, = ac/y is a ‘Reynolds’ number 
_ based on the speed of light. As an example pooeide? the case of mercury where 
A= 10-7, M=100 (the experimental upper limit) and R, = O(10") (with 
aa = 1 ten’ In this case (43) becomes V, ,,/¢c < 10’ which is obviously always 
satisfied. 

4 It should not be concluded from the previous paragraph that the effect of 
@, 1s completely negligible, for it has a dominant effect on the electric field. The 
free charge in the boundary causes the electric field to make a large increase 
on passing through this layer. This is best illustrated by a concrete example. 
Consider a long rectangular duct with non-conducting walls. Except in the 
_ boundary layer near the ends the solution is essentially that given by equations 
(14) and (15), with gy = 0. Using this solution we can calculate E from equation 


Le= 3 Ho ror + V2 Bo, 4) 
1 OB, 
Pee a, a 


LA a Op | 1 i 
eae mae (45) 
p= 0 


It is seen that E, is constant. One might suppose that F, is constant all the 
way to the end of the duct, taking a jump through the surface of the wall. This 
is not the case, since for nonconducting walls B, is zero at the wall, therefore 
0B./Oy is zero on the vertical wall, and equation (44) then shows that EF, = 0 
on the vertical wall. We conclude that E, varies from zero at the wall to the 
value given by (45) as we go through a boundary layer of thickness alVM : 
This is due to the free charge in this boundary layer. (There is some question 
as to whether the free charge causes the variation in electric field or whether 
the variation in electric field causes the free charge.) 
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Conclusion 


In the flow of conducting fluids through ducts, in general, the domain of 
the equations describing the flow is not the same as the domain of the fluid. It 
was found that the problem reduced to two sets of equations, one set, (11) and 
(12), in the fluid, and one set, (13), in the wall, with boundary conditions 
specified on the outer boundary of the wall and across the inner boundary 
of the wall. The known special cases (SHERCLIFF [3, 4]) where this problem 
reduces to an ordinary boundary value problem with boundary conditions given 
at the boundary of the fluid, follow as limiting cases of the above formu- 
lations. These are y > 0, » > oo andh > 0. 

In the bulk of the paper the flow through rectangular ducts was considered. 
The essential conclusion is that increasing the wall conductivity tends to 
decrease the average velocity if the pressure gradient is the same in the two 
cases. Stated differently, when the wall conductivity is increased, the pressure 
gradient must be increased in order to maintain the same mass flow. 

The flow through arbitrarily shaped symmetrical ducts was considered for 
large M and arbitrary wall conductivity. It was found that when » > oo the 
velocity tends to be uniform except in a thin boundary layer along the wall. 

In the final section the effect of the neglected free charge was considered. 
We found that free charge tends to accumulate near walls which are parallel to 
the applied field, but not enough to effect the velocity distribution. 
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Zusammenfassung 


In dieser Arbeit wird die Strémung einer elektrisch-leitenden Fliissigkeit durch 
ein gerades Rohr mit einem gleichférmigen, querlaufenden magnetischen Feld 
betrachtet. Das Problem wird unter Beriicksichtigung des elektromagnetischen 
Feldes sowohl innerhalb als ausserhalb der Fliissigkeit formuliert. Besondere Auf- 
merksamkeit wird auf die Ableitung der Randbedingungen gerichtet. Es wird 
klargemacht, dass, wenn die Wande des Rohres endliche Leitfahigkeit haben, das 
Problem, abgesehen von Einzelfallen, kein gewohnliches Randwertproblem dar- 
stellt. Eine Losung wird fiir die Strémung durch ein rechteckiges Rohr mit unend- 
lich fernen, dem ausseren Feld parallelen Wanden bei beliebiger Wandleitfahigkeit 
gefunden, ebenso fiir die Str6mung durch ein endliches rechteckiges Rohr mit ideal- 
leitenden Wanden. Es gelingt, die asymptotische Geschwindigkeitsverteilung bei 
grosser Hartmannscher Zahl fiir ein diinnwandiges symmetrisches Rohr bei belie- 
biger Wandleitfahigkeit anzugeben. 
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The Influence of Turbulence on the Transfer of Heat to 
Cylinders near the Stagnation Point 
By J. Kerstin"), P. F. Marper!), and H. H. Socrn?), Providence, R. I., USA 


1. Introduction 


In a previous paper [1]*) two of the authors have shown that the intensity 
_ of turbulence of a free stream exerts an important influence on the rate of heat 
_ transfer from a cylinder in cross-flow. The most noteworthy conclusion from 
_ that investigation was the realization that the effect was a local one and 
_ additional to the expected influence from the effect of free-stream turbulence 
on transition and separation. In another paper [2] it was shown that this local 
_ effect was absent in the case of a flat plate at zero incidence, but that it re-ap- 
_ peared when a favorable pressure gradient was imposed on it. In both cases, 
_ the local Nusselt number increased by amounts of the order of 25-50°/4 when 
_ the turbulence intensity increased relatively little, say from 0-5 to 2%. A 
_ plausible explanation for these effects was given elsewhere [3], it being conjec- 
- tured that the local Nusselt number increased as a result of the oscillations 
carried by the free stream; it was, namely, shown that any oscillations accom- 
panied by a pressure gradient must be expected to modify the mean velocity 
profile of an associated boundary layer. Changes in the velocity profile will, in 
turn, modify the temperature profile and hence the Jocal Nusselt number. 
Although the preceding explanation appears plausible, the details of the 
mechanism operating under such conditions are far from being clarified and 
understood. Nevertheless, if the explanation is plausible, one must expect 
large local increases in the rate of heat transfer in the neighborhood of the 
forward stagnation point of a cylinder (or any bluff body) in cross-flow when 
the boundary layer is laminar and where the free-stream velocity gradient is 
large, and where, therefore, there exists a large downstream gradient in the 
amplitude of oscillation, as required by the theory advanced in reference [3]. 
The present paper addresses itself to the experimental verification of this point. 
The increase in the local rate of heat transfer from a bluff body which is 
produced by an increase in the free-stream turbulence intensity has been 
observed before. Even more, it appears that evidence of this effect can be 


1) Department of Engineering, Brown University. 
2) Department of Mechanical Engineering, Tulane University, formerly with Brown University. 


3) Numbers in brackets denote the References, page 131. 
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discerned in experimental investigations which were carried out at a time 
when the importance of the turbulence in the free stream in general had not 
yet been discovered. ae 

Figure 1 shows a comparison between the classical experiments performed 
by Scumipt and WENNER [4] and the theoretical calculations due to Fross- 
LING [5]. Since the resulting variation of the Nusselt number with the angular 
co-ordinate @ is fairly sensitive to the assumed variation of free-stream velocity 
U with the co-ordinate x = Rq@ (y= radians), the equation given by FROss- 
LING was re-evaluated with the aid of the following series expansion for U, 
where U,, is the velocity of approach, D being the diameter of the cylinder: 


Pe Ue {3-631 (4) - 3.275 (4) — 0-168 (5)}- (1) 


The coefficients in equation (1) were given by ECKERT [6] who deduced them 
from SCHMIDT and WENNER’S measurementsat a Reynolds number Re = 170000. 
Thus we were led to the expression 


Nu 


x \2 x \4 
Tare = 09449 — 0-7693 (+) — 09-3009 (=)  (Pr=0-7) 2) 


in which the Reynolds and Nusselt numbers are referred to the diameter D. The 
group Nu/ Re? is very convenient for presenting experimental results, because 
it is independent of Reynolds number and is a unique function of the angular 
coordinate for a given Prandtl number, Py, as long as the coefficients in the 
expansion in equation (1) have been chosen properly. 

The diagram presents three series of results, each for a different diameter D 
but for similar ranges of velocities U,,. Hence the ranges of Reynolds numbers 
are also different. For the sake of simplicity, the same theoretical curve, equa- 
tion (2) is shown in all three graphs, but strictly speaking, the power-expansion 
in equation (1) most closely corresponds to the diagram in Figure 1c. In any 
case its accuracy becomes questionable at g = 60°. 

From the nature of the theoretical calculation which neglects free-stream 
oscillations, it is clear that it refers to a stream with zero intensity of turbulence, 
i. e. to a laminar stream. The diagrams show that all measured points indicate 
a higher rate of heat transfer for the range of the angular coordinate 
(0 < m < 60°) where the boundary layer is expected to be laminar, as compared 
with the theoretical curve. The fact that all points lie above the theoretical 
curve suggests that the deviations are not due to random experimental errors. 
Furthermore, at the highest Reynolds numbers (Re = 52800, 170000, 426000 
respectively) the Nusselt numbers at the stagnation point exceed the calculated 
value by 7%, 13% and 16% approximately which is considerably larger than 
the 5% margin of error of the determination. This suggests that the deviations 
are systematic and increase as the Reynolds number is increased. In actual 
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Local rate of heat transfer from forward portion of a cylinder. Comparison between experiments 
due to Scumipt and WENNER [4] and calculations due to FrRoEssLinG [5], Pr = 0:7, free-stream 
velocity from equation (1). 

a 8290 < Re< 52800, Di= 5 em. 

b 15550 < Re < 170000, D = 10cm. 

c 39200 < Re < 426000, 1 D=25 em. 
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fact, as was shown in reference [1], it is highly probable that the free stream _ 
available to SCHMIDT and WENNER was one of relatively low turbulence — 
intensity (of order 0-9%), but also one in which the intensity of turbulence 
increased with the velocity, and hence with the Reynolds number. Thus the 
apparent increase in Nusselt number with the Reynolds number can possibly — 
by interpreted as an increase with increasing intensity of turbulence. : 

If this hypothesis is correct, then the curves for a given velocity, and hence 
for a given intensity of turbulence, should most nearly coincide, the only 
difference being now due to the different pressure distributions, as expressed 
by U(x) and associated with different cylinder diameters. 

The results for the two groups of Reynolds numbers, 52800 in Figure 1a, 
101300 in Figure 16, and 257600 in Figure Ic satisfy this condition reasonably 
closely. The two sets of results have been replotted in Figure 2 whose appear- 
ance confirms our guess, in so far as the general spread in the values has now 
decreased. 

Figure 3 shows a comparison between the measurements performed by 
five groups of investigators [4] and [7-10] all at about the same Reynolds 
number of Re = 39600. For illustrative purposes, FROSSLING’sS curve from 
equation (2) has also been drawn. Once again all measured values exceed the 
calculated ones and the excess at the stagnation point reaches an order of 
magnitude of 60-70%. It is unlikely that so many experimental investigations 
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Verification for two sets of three Reynolds numbers which correspond to approximately equal 
velocities (and turbulence intensities). Replot from Figure 1. 

a Re=52800, D=5cm; Re = 101300, D=10cm; Re = 257600, D = 25cm; 
b Re = 21200, D=5cm; Re= 39800, D=10cm; Re = 102000, D = 25 as 
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could show such large errors, all of them positive, and it is probable that they 
are due to systematic differences in the recorded intensities of turbulence in 
the various streams. 
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Figure 3 
Variation of local Nusselt number with angular co-ordinate for cylinder in cross-flow at Re = 39600 


_ (or a value close to it). Measurements by five groups of investigators, theory by FrOss.Linc [5]. 


a Kevyn [7], Re = 39500, d Kren [10], Re = 39 600, 
b Sma tt [8], Re = 39600, é ScumiptT and WENNER [4], Re = 39800. 
c Drew and Ryan [9], Re = 39600, 


2. Scope of the Investigation 


During the present experimental investigation, measurements were made 
of the local coefficient of heat transfer from a vertical cylinder placed in a wind 
tunnel. The measurements were made at varying velocities and with intensities 
of turbulence which were increased by placing suitable screens upstream of the 
model. Hence the two parameters, the Reynolds number Re, and the intensity 
of turbulence, Tu, were varied nearly independently. A completely independent 
adjustment proved to be impossible because the intensity of turbulence pro- 
duced by a given screen turned out to be a function of air speed [1]. Finally, 
measurements were made over the full range of angular positions g, from 0° 
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to 180°, but owing to the complexity of the results, and the difficulties in their 
interpretation, the present paper will discuss exclusively the case of the 
laminar boundary layer, say up to y = 60°. 
The measurements covered a range of Reynolds numbers 
Us. D 
v 


120 x 10°?< Re= < 290 x 10? 


and a range of turbulence intensities (based on the longitudinal fluctuation w’) 


(u’2)1/2 


05% <e Tu = at, Sd a ZL 2 


3. Description of Experimental Arrangement 


The model used in the present investigation consisted of a 4-21” OD brass 
tube provided with a 1/2” x 10” strip heater, similar to that used by SCHMIDT 
and WENNER [4]. The tube was placed vertically inthe 20” x 32” Brown 
University low-speed tunnel, [1], and could be rotated about its vertical axis. 
The angle of turn was set with the aid of two indexing plates accommodated 
outside the tunnel test-section. 

The design of the tube is depicted in Figure 4 which contains two cross- 
sections and a view of the model. The heat transfer tube B was placed in a slot 
in the brass tube A, and consisted of a specially wound, sheathed nichrome 
spiral C cast into a copper slab. The copper slab was enclosed in a housing D 
milled out of a solid brass bar. The housing insulated the element from the 
interior of tube A, and a brass pipe F provided an outlet for the heater leads, 
as well as for several thermocouples. A micarta slab F connected the heat 
transfer element to the housing, insulating it thermally. The thin gap C be- 
tween the tube wall and heater slab was filled with a plastic, heat resistant paint. 
The housing D together with the element B were joined to tube A with screws, 
and the whole assembly was turned to size and polished in a lathe. Upon 
heating, the heater element had a tendency to protrude somewhat, but the 
resulting unevenness was less than 0-001”. 

The interior of tube A was supplied with saturated steam from a small 
electrical boiler, and the heater C was supplied with DC current of a highly 
stabilized voltage. During an experiment, the current supplied to the heat 
transfer element was so adjusted as to keep its surface temperature substan- 
tially equal to the surface temperature of the tube. At low speeds the maximum 
difference in temperature did not exceed 0:5°C rising to 2°C at the highest 
wind speeds. In this manner, nearly all the Joule heat evolved in the electric 
heater was forced to leave through the exposed surface of the heat transfer 
element, and was inhibited from leaking out sideways by the small temperature 
difference coupled to a high thermal resistance. Sample calculations showed 
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that possible leaks did not exceed 5% of the quantity measured. Since no 

greater accuracy was intended, no correction for heat leaks was applied to the 

final result. Thus, in contrast with most other investigations [7, 11, 12, 13], 
. ‘the measurement of the rate of heat transfer was direct. Strictly speaking, the 
“Measurement extended over a width of 1/2” (14-5° of angle), and thus contri- 
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The experimental tube. 


‘buted a mean measurement over a fraction of the circumference. It can be 
easily shown that for the measured variation of Nusselt number with angular 
coordinate (convex upwards), this averaging tended to underestimate the rate 
of heat transfer. 

The electrical input was measured with the aid of a sensitive potentiometer, 
as described in reference [1]. The temperatures were measured with the aid of 
copper-constantan surface thermocouples whose location is shown in Figure 4. 
The wires were led out of the assembly through a tube shown in the drawing. 
In addition, provision was made by two radiation-shielded thermocouples, one 
inside cylinder A to measure the steam temperature, and one in the settling 
chamber of the wind tunnel to measure the air temperature. All temperatures 
were measured differentially with respect to the steam junction, and the EM F’s 
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were recorded on a Brown potentiometer-recorder using a scale of 1 mV to 
25 cm. In all, temperatures were measured with an accuracy of 0-15°C. 

The tube was provided with two static-pressure probes H, Figure 4, placed 
symmetrically with respect to the center plane of the heater; they were used to 
position the tube with respect to the wind. 

All remaining details and procedures, including those for the evaluation of 
the results were identical with the ones described in reference [1] and need 
not be repeated here, except for re-stating that the Nusselt and Reynolds 
numbers were evaluated with thermal properties averaged integrally over the 
boundary layer thickness, as in reference [4]. 


4. Experimental Results 


The experimental results cover measurements in three distinct, non-over- 
lapping ranges of turbulence intensity. The ranges were obtained in the free 
tunnel and by putting two, different, turbulence-producing screens, 12 inches 
upstream of the model. As already mentioned, neither in the free tunnel, nor 
in the presence of the screens, did the turbulence intensity remain constant 
over the range of wind velocities available in the tunnel and it was necessary to 
calibrate the turbulence intensity in terms of air velocity. The calibration 
curves, reproduced from reference [1], are shown in Figure 5 for easy access. 
The resulting ranges of parameters covered in the three runs are shown summa- 
rized in Table 1. It will be noticed that the runs were arranged to avoid over-. 
lapping in turbulence intensity so that at any Reynolds number a point in 
run 1 had the lowest, a point in run 2 had an intermediate, and a point in run 3 
had the highest intensity of turbulence. 


Table 1 
Ranges of parameters covered 
Secen Range of 
Gariee Desig: Wire titans Range of Paya numbers, 
nation| Mesh | dia- intensity, Tu é 

meter 
Run 1 | 
Low Intensity . . . . | None | None|None| 1-0 to 0-5 | 132 x 108 to 290 x 103 
Run 2 
Intermediate Intensity | No. 2} 4/,” |0:062”| 1:9 to 1-5 130 x 10® to 260 x 103} 
Run 3 
High Intensity. . . .|No.1| 3/," |0.148”| 2-1 to 2-7 | 120 x 103 to 220 x 108 


tn carrying out the experimental program, it was found impracticable to 
adjust the tunnel to a given speed and to make measurements over the whole 
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range of angular positions. Instead, the heater was turned to a given position, 
_and measurements were performed over a range of Reynolds numbers. In all 
_ cases, measurements were made at two symmetrical positions, denoted by + 
-and — g respectively, and averaged graphically. Owing to the tunnel charac- 
_ teristic, a certain lack of symmetry was unavoidable. 
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Figure 5 


Turbulence intensity in terms of air velocity from reference [1]. 
a) Low intensity, no screen, 
b) Intermediate intensity, screen No. 2 at 12”, 
c) High intensity, screen No, 1 at 12”. 


The direct experimental results are seen plotted in Figure 6. They are 
plotted in the form of a diagram of the group Nw/Re'? because for a given 
Prandtl number, this ratio is a unique function of the angular coordinate at 
Tu = 0, as already stated earlier. It is, however, necessary to remember that 
equation (2) remains valid only for a given pressure distribution near the 
stagnation line. Hence the variation of Nu/Ret/? with Reynolds number may 
be due to two causes. First, as the Reynolds number varies, the pressure 
distribution may vary somewhat, causing a variation in Nu/Re1!; it is supposed 
that this is a small variation. Secondly, as the Reynolds number varies, the 
turbulence intensity varies, and affects the value of Nu/Re1!; it is supposed 
that this constitutes the major effect. 
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Figure 6 


Experimental results plotted as the ratio Nu/Re!!? versus the Reynolds number Re, at different 
angular positions + @. 
a Run 1, lowest turbulence intensity, 
b Run 2, intermediate turbulence intensity, 
e Run 38, highest turbulence intensity. 
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An idea about the effect of different pressure distributions can be formed 


with reference to Figure 7 which contains plots of the variation of the heat 


_ transfer parameter Nu/Re1!? with the angular coordinate m for two pressure 
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Figure 7 


Effect of pressure distribution on heat transfer parameter Nu/Re!/? 
a Re = 170000, equation (2), b Re = 19000, equation (3). 


distributions. The lower curve is a plot of equation (2) which corresponds to 
the pressure distribution implied in equation (1) and obtained experimentally 


at a Reynolds number of Re =170000. The upper curve corresponds to 


FROSSLING’S [5] explicit equation 


ae = {0- 9449 — 0-5100 (4.)° — 0:5956 (5). (Pr=0-7) (3) 


obtained from the pressure distribution measured by HIEMENz [4] at a Rey- 


_nolds number Re = 19000 and corresponding to the velocity distribution 


Ulett fs, {3-6314 (4) — 21709 (+ -) — 15144 (5)}- (4) 


It is seen that the two results agree very closely near the stagnation point, 
diverging as the angular coordinate m increases. The discrepancy reaches a 
value of 6:5°% at m = 60° which means that the preceding argument ceases to 
be valid there. 

The two Reynolds numbers of 19000 and 170000 extend over a sufficiently 
large range to make it certain that in the range of Reynolds numbers covered 
in the present investigation, the possible discrepancies from this source would 
be smaller than those suggested by the diagram in Figure 7. 
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In order to exhibit the local effect of turbulence intensity, it is best to 
proceed by cross-plotting the results from Figure 6 for several Reynolds num- | 
bers. At a constant Reynolds number, and for a given screen configuration, the — 
intensity of turbulence is constant. Since the difference between the three runs 
consists only in the insertion of a screen ahead of the model for runs 2 and 3, the 
pressure distribution will remain virtually unaffected near the stagnation region, - 
as confirmed by a subsequent investigation [14]. If there were no local influence — 
from turbulence intensity, all such cross-plots should arrange themselves around 
the theoretical curve, equation (2), in a band corresponding to the accuracy of 
the measurement. In the contrary case, they should form a distinct family of 
curves, with the intensity of turbulence as a parameter. Three such crossplots 
are shown in Figure 8. The intensity of turbulence corresponding to each point 
plotted was obtained with reference to the calibration curves given in Figure 5. 
The theoretical curve from equation (2) has been drawn in for comparison; it 
corresponds to Tu = 0. 

The diagrams unmistakably demonstrate the local effect of turbulence 
intensity on the rate of heat transfer and hence on the characteristics of the 
associated laminar boundary layer. The effect is unexpectedly large, reaching 
a value of 80% at the stagnation point for a change of turbulence intensity 
from 0% (calculated) to less than 3% (measured). 

The same results are seen plotted in an alternative way in Figure 9 in which 
the group Nu/Re1/2 has been plotted against Reynolds number for different 
values of the angular co-ordinate g. The important feature to be noticed in this | 
diagram is the fact that, once more, the experimental points do not group 
themselves around the theoretical constant values shown in broken lines. It 
will be remembered that along each curve the turbulence intensity, and pos- 
sibly the pressure distribution, vary somewhat with Reynolds number. How- 
ever, owing to the fact that the three ranges of turbulence intensity do not 
overlap, the local effect on the rate of heat transfer exerted by turbulence 
intensity all along the laminar boundary layer is apparent. At one given 
Reynolds number, particularly near the stagnation point, the pressure distri- 
bution is insensitive to turbulence intensity, and the diagram clearly demon- 
strates that the heat transfer parameter N u/Re!? increases from run to run, 
as the turbulence intensity is increased. 

In order to gain an insight into the effect of turbulence intensity, it is 
convenient to cross-plot the experimental results from Figure 6 in yet another 
form. The diagram in Figure 10 shows the variation of the ratio N u/Re ata 
given intensity of turbulence to that at zero intensity, denoted (Nu/ Retl?),, at 
different angular co-ordinates y and for different Reynolds numbers. It is seen 
that the effect is larger at lower intensities of turbulence and nearer the 
stagnation point and that it depends to some extent on the Reynolds number. 
Hence it may be concluded that a given increase in the turbulence intensity 
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Figure 8 
Variation of local Nusselt number on a circular cylinder with turbulence intensity Tu and angular 
co-ordinate m. Crossplot from Figure 6. 
@ Re= 140 x 103; b Re'= 180° x 105; 6 Re = 220 x 10°; 
----- theory after FRoEssLING [5], equation (2). 
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Variation of local Nusselt number on a circular cylinder with turbulence intensity Tu and angular 
co-ordinate gm, Crossplot from Figure 8. Runs 1, 2,3 correspond to non-overlapping bands of 
increasing turbulence intensity. 

a p=0, b pM = 20°, c p=40°, d p=60°, 

—----- theoretical value after FRoEssLiNG [5], equation (2). 
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Variation of the ratio (Nu/Re!2)/(Nu/Re1!?),) of the heat transfer parameter at given intensity of 
turbulence to that at zero intensity of turbulence with turbulence intensity Tu at different angular 
coordinates m and different Reynolds numbers Re. (Nu/Ret!?), denotes the theoretical value of the 
ratio at given angular coordinate m from equation (2) i.e. for Tu = 0. 
a p=, b p= 20°, c p=40°, d p= 60°. 
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produces a larger increase in the Nusselt number in regions where the gradient 
of free-stream velocity is larger. 


5. Concluding Remarks 


The preceding, preliminary investigation demonstrates clearly the existence 
of a local effect from turbulence intensity on the Nusselt number across a 
laminar boundary layer with a favorable pressure gradient imposed on it. The 
effect appears to be largest at low intensity of turbulence and at points asso- 
ciated with the largest gradients of free-stream velocity. 

The present measurements are not considered to have been of a sufficiently 
high precision in order to permit us to supply definite, quantitative measures 
for this effect, but the magnitude of the effect was sufficiently large to justify 
the preceding qualitative conclusions. 

The preceding discussion, admittedly, constitutes a somewhat simplified 
interpretation of the experimental results. This is due to the fact that the 
changes in the structure of the free stream were not investigated in detail in the 
distance between the entrance section and the model. It is well-known that the 
oscillations carried by a free stream are amplified along streamlines, the largest 
amplification occurring along the stagnation streamline. It is, therefore, clear 
that the modification in the velocity profiles in the laminar boundary layer 
caused by oscillations will be governed by their amplitude outside the boundary 
layer, and not by those in the free tunnel. However, the detailed mechanism 
of amplification is not known, and one can only surmise that the increase in 
the heat transfer parameter N. u/Re? with Reynolds number evident in 
Figures 6 and 8, even for the clear tunnel and for screen No. 2 for which the 
intensity of turbulence decreases with Reynolds number, is a consequence of 
this phenomenon. It is clear that a complete explanation cannot be obtained 
other than by further detailed measurements. 

It is pertinent to remark here that in a subsequent investigation on the 
effect of turbulence intensity on mass transfer from a cylinder in cross-flow, by 
H. H. Socin, V. SUBRAMANIAN and R. J. Socin [14] confirmed the preceding 
qualitative results. Quantitatively they found that the rate of mass transfer 
at a turbulence intensity of the order of Tu — 1-0% agreed quite well with the 
theoretical values showing a much smaller effect in the range 0-1%. By a suit- 
able choice of average values of the transport properties, namely by forming 
the Reynolds and Nusselt numbers with values of viscosity and conductivity 
corrrsponding to the mean boundary layer temperature, and by inserting the 
free-stream density into the Reynolds number, it is possible to bring the 
values of the heat transfer parameter Nu/Re? for the clear tunnel to good 
agreement with equation (2). The essential point, however, is that a further. 
increase in turbulence intensity causes large further increases in the Nusselt 


—_ eta 


Vol. XII, 1961 The Influence of Turbulence on the Transfer of Heat to Cylinders 131 


“number. Hence it may be concluded that the particular manner in which the 
Properties Lt, k and @ are averaged, has no bearing on the qualitative conclusions 
advanced in this paper. 
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Zusammenfassung 


Es wird das Problem des Einflusses der Turbulenz auf den Warmeiibergang an 
einem quer angestromten Zylinder untersucht. Die Untersuchung beschrankt sich 
auf den Teil des Umfangs, iiber den sich die laminare Grenzschicht erstreckt. Einige 
altere Messungen an der Vorderseite eines Zylinders werden miteinander verglichen, 
und es zeigt sich, dass ein wichtiger Parameter, namlich die Intensitat der Vortur- 
bulenz bei den Untersuchungen nicht beriicksichtigt worden ist. Die bekannten 
Messungen von Schmidt und Wenner werden kritisch untersucht und der genannte 
Effekt wird an ihnen aufgezeigt. 

Weiter werden eigene Messungen der lokalen Nusseltschen Zahl beschrieben. Die 
Messungen waren durchgefiihrt worden an der Vorderseite eines einzelnen Zylinders 
im Querstrom, und zwar in drei sich nicht iiberdeckenden Gebieten der Intensitat 
der Vorturbulenz. Die Versuchsanordnung wird beschrieben und mégliche Einfliisse 
der Veranderung der Druckverteilung am Zylinder werden diskutiert. Obwohl dieser 
letztere Effekt nicht véllig eliminiert worden war, zeigen die Versuche doch klar, 
dass steigende Turbulenz ein starkes Anwachsen der lokalen Nusselt-Zahl bedingt. 
Dieses Anwachsen ist bei kleiner Turbulenz besonders stark. Ein Vergleich mit der 
theoretischen Losung Froesslings, die fiir Turbulenz 0 gilt, zeigt, dass die Turbulenz 
den Warmeiibergang in der Staulinie bis zu 80° erhéhen kann, bei einer Turbulenz- 
variation von 0-3%. 

Diese Ergebnisse, die qualitativ interpretiert werden, stiitzen die Hypothese 
einer friiheren Veréffentlichung, die annimmt, dass der Effekt durch Anderungen 
der laminaren Grenzschicht, die hervorgerufen werden durch Schwankungen von 
zunehmender Amplitude, wie sie bei zunehmenden Vorturbulenzen auftreten, 
bedingt ist. 


(Received: September 22, 1960.) 


Biharmonic Solutions to the Steady-State Thermoelastic 
Problems in Three Dimensions | 
By Jerzy Nowrnskt, Madison, Wis., USA}) 


(An abstract of this paper has been presented on the International Con- 
ference on Partial Differential Equations and Continuum Mechanics in Madi- 
son, June 1960). 


1. Introduction 


Because of a rising interest in thermal problems, much recent literature has 
dealt with the complex question of thermoelastic fields in three dimensions. 
Thus, confining ourselves to the case of a steady state of temperature in an 


1) Army Mathematics Research Center, University of Wisconsin. 
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isotropic body, we may cite the feoresent tivehiork by MINDLIN andCHENG [1]?), 
TERNBERG and McDowELt [2], SNEDDON and Locket [3], and McDowELL 
and STERNBERG [4]. The first three papers treat the semi-infinite body?), while 
the last concerns a thick spherical shell in an axisymmetric temperature field‘). 
The diversity of methods used in these papers, even for related problems‘), 
_ suggests that it would be expedient to look for still another method having a 
_ somewhat broader applicability. One such method already has been introduced 
a in ordinary (isothermal) three-dimensional elasticity and is due apparently to 
_ KELvin [8]®), and ArmAns1 [10]. It has been used extensively by TREFFTz [11]. 
_ The method corresponds, to some extent, to the method of KoLosov-MusKHE- 
_LISHVILI, based on the Goursat representation of plane biharmonic functions, 
A Behich has proved to be so effective in deducing explicit solutions to many 
4 _ important problems in two dimensions. 
c In what follows we deal with the application of the Kelvin-Almansi method 
- to thermoelastic problems in three dimensions concerning infinite, semi-infinite 
_ and spherical regions. To make the paper self-contained we first briefly explain 
_ the basic equations, for future reference. 


a 


Do 


AERP PRE 


a a vee 


a 


2. Basic Equations 


Consider an isotropic homogeneous elastic medium occupying a region of 

_ space D with the boundary B. Let the body be in static equilibrium subject to 
| a steady-state temperature field 7(P), P being a point in D + B and T being 
measured above a uniform temperature in which the body remains unstrained. 
With no loss in generality, we may assume that surface tractions and body 
forces are absent, since by virtue of the linearity of the problem, postulated in 
the following, the influence of an external load can be determined by solving 


an ordinary elastic problem’). 
Within the classical theory of elasticity, the pure thermoelastic problem 
(in the aforementioned sense) is governed by the linear strain-displacement 


2) Numbers in brackets refer to References, page 148. 

3) In the paper [3] the case of an infinite thick plate is also solved. Apparently, this case was 
previously solved by Lurie [6], page 193, using thermoelastic potential, and SHarma [7] using 
Fourier integral representation. 

4) The limiting cases being a solid sphere and an infinite medium with a spherical cavity. 

5) In [1], an arbitrary distribution of temperature is treated by an extended scheme of Goo- 
DIER [5]; in [2] a steady-state heat flow is approached by the method of Green; and in [3] double 
Fourier transforms are used. In the paper [4] the method of thermoelastic displacement potential 
is employed. 

6) See in this connection the book by SoxornikorFF [9], p. 335 and 351. 

7) It is implied that the temperature increase T(P) does not effect the value of the elastic 
moduli and the coefficient of thermal expansion. In the last case the third term in the Duhamel- 


5 ry 
Neumann equations of state (2) takes the form | u ae | «(9) dd. 
; ft) 


| 
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and stress-strain relations together with the homogeneous equations of equi- 
librium. With reference to a rectangular cartesian coordinate system *; 
(i = 1, 2, 3) and in index notation, these conditions provide the following field 
equations: 


v 1+y 
t= 24 [tug + (Goay ee Toy 7) Bs, 2) 


where the usual conventions of summation (by repeating indices) and differen- 
tiation (by means of a comma) have been used. ;, e;; and t;; denote, as usual, 
the components of displacement, strain and stress, respectively; 6;; is the 
Kronecker delta, w the shear modulus, y Porsson’s ratio and « the coefficient 
of thermal expansion. 

To the general field equations (1), (2), (3), which are valid throughout the 
region D occupied by the body, we must add the boundary conditions which, 
for the pure thermoelastic problem (i. e., for the boundary B free of tractions) 
take the form Bey Ty er (4) 


UP ae | 


vy; being the components of a unit vector normal to the surface B. 

Substitution of (1) and (2) into (3) provides the Duhamel-Neumann field 
equations in terms of displacement. 

After changing from indicial notation to vector notation, which is more 
expedient at this stage, we obtain 


Au + x grad divu — B grad T=0, (5) 


where wu denotes the displacement vector, E YouNG’s modulus, A the Laplace 
space operator and x= 1/(1— 27), B=Eaxju=2 (1+) a/(1—2>). 

Let us now turn to the temperature field. Assume that the heat is generated 
at no more than a finite number of points situated in D+ B. Then for a 
steady flow of heat, the governing temperature field equation is the Laplace 


equation 
ATO, (6) 


provided that the body is thermally isotropic and homogeneous. Thus T(P) is 
a harmonic function, i. e. regular throughout the region D + B, except at the 
points of generation of heat where it may become singular§), 


Denote the dilatation by 
G=divu (7) 


8) By the boundary value or the value of a function at a singular point we understand in the 
sequel the value of the function obtained by a suitable limit process. 


PET aS 
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_and take the divergence of both members of (5). Because of the temperature 
_ equation (6), this gives 


. 


»? 


Ab =0, (8) 


: a relation which is well known for the isothermal problem. Now, taking the 
Laplacian of both members of (5), we derive, in view of (8) 


4 3 AAé = 0. (9) 


_ Thus every component of the displacement field, in the steady-state thermo- 
elastic problem under investigation, represents a biharmonic function, as in 
_ the isothermal case. This permits us to utilize in the solution of the thermo- 
_ elastic problem a method which was worked out for the isothermal case. 

In preparation for an application of the method) let us observe that if 
_¢(P) and y(P) are two arbitrary functions harmonic inside the region D and 
if y is the modulus of a position vector 7 = x, %9, x3, then 


P(P) = o(P) +? y(P) (10) 


_is a solution of the biharmonic equation in the region. This assertion may be 
proved by direct substitution of (10) into (9). Moreover, the converse is true?®). 

_ Every biharmonic function Y/(P) in D can be expressed in the form (10), where 
¢(P) and w(P) are harmonic in D. Thus we have a general representation of 
biharmonic functions in terms of pairs of harmonic functions. Since Equation (9) 
was derived from Equation (5) by differentiation we can assert that every 
solution of (5) is a biharmonic function, but the converse may not be true, i. e., 
not every biharmonic function is a solution of (5). In order to obtain a represen- 
tation formula for the solution of the differential equation (5) we must fulfill 
the additional requirement (8), viz. that the divergence of the biharmonic 
vector function “ must be harmonic function. 

Represent the displacement vector in the form 


=> 


u(P) = 6(P) + 7 grad y(P) , (11) 


where 4(P) is a harmonic vector (a vector whose components are harmonic 
functions) and y(P) is a harmonic scalar function. Since the components of the 
gradient of a harmonic function are themselves harmonic, the representation 
(11) provides, by virtue of what was said in connection with (10), s solution of the 
field equations (9). It may be shown by simple calculation that (11) also becomes 
a solution of (5) if the following relation between the harmonic functions »p(P) 


and $(P) exists: iia he 
y 2+ % CET 1 
v= oe D(P) dr. (12) 
0 


§) See ALMANSI [10], TreFrrz [11], BERGMAN and SCHIFFER [12]. 
10) See [12], page 230. 
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Here the integration is to be performed along a fixed position vector and 
0 hegre 
DP) =2[y+ 2+ arze|=—xdivg+ PT +C. (13) 


The constant of integration C which depends on the direction of integration 
has to be determined in such a way that D(P) really is a harmonic function”). 
Then, by a known lemma?2), y will also be harmonic, which implies the har- 
monicity of the divergence of # in (11), as was demanded. It is now possible 
to construct an infinite set of particular solutions of (5) by starting with an 


arbitrary harmonic vector function d (i.e. with three arbitrary harmonic 
scalar functions ¢;, 7 = 1, 2, 3) and deriving the fourth harmonic function y 
by means of the relation (12). 

We now observe that since the stresses arise from the displacements by 
simple differentiation and the operator of partial differentiation commutes 
with the biharmonic operator, the stresses induced by the steady heat flow 
are biharmonic functions, a result known in ordinary elasticity. 


3. Some Theorems on Harmonic Functions 


As a further preliminary we shall quote, for future reference, two theorems 
concerning harmonic functions. The regions that we shall consider here are the 
half-space, the plate and completely regular regions. The latter are defined as 
closed finite or infinite regions that are bounded by closed smooth surface?%). 
In the case of infinite regions, we shall assume that the harmonic functions 
considered are regular at the point at infinity. This means that the functions 
tend to zero at infinity no slower than 1/o and their partial derivatives no 
slower that 1/07, where @ designates the distance from any fixed point. 

Bearing this in mind we formulate the following theorems. 


Theorem 1. A function which is harmonic in a completely regular region and 
vanishes at all points of the boundary must vanish identically in the region. 


Theorem 2. A function which is harmonic in a half-space and vanishes at all 
points of the bounding plane must vanish identically in the half-space. The proof of 
the first theorem may be found in any book on potential theory!4). A proof of 
the second theorem is given in the Appendix. 


11) Of course the middle expression in (13) is harmonic, since yw is harmonic by hypothesis 
and from this results the harmonicity of r (Ow/0r) as may be trivially proved. Clearly, the right-hand 
side of (13) is harmonic if C is harmonic, since the divergence of a harmonic function is harmonic 
and T is harmonic by (6). 

12) [12], page 230. 

18) A smooth surface is a surface with continuously turning tangent plane. An example of an 
infinite completely regular region provides the entire space exterior to a spherical surface. Our 
requirement of complete regularity of the region may be considerably weakened. See in this con- 
nection the regularity conditions introduced, e. g., in [13], page 113 and 217. 

14) ‘See e. g. [13], p. 213. 
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We now proceed to ie our general considerations to the steady-state 
thermoelastic problem for some fundamental regions. 


4. The Half-Space 


___ Denote now the cartesian coordinates by x, y, z. Let the half-space occupy 

the region z > O with the boundary z = 0. 

__ We shall show that for a steady flow of heat, s simple particular solution 

of the type (11) permits us to obtain a general solution for the case considered. 
To this end drop the nonessential constant in (13) and assume in addition 

that = 

divg=2T. (14) 


Then from (13) we find that D(P) =0 and from (12) that y(P) =0. This 
shows that we seek the solution to the problem in the form 


U,v,W = #1, bs, bs , (15) 


where u,v, w denote the components of the displacement vector uw and ¢; 
(4 = 1, 2, 3) are harmonic components of d. Since from (7) in the case con- 
sidered 

0=divd, (16) 


we confirm the fulfillment of the condition (8) that the dilatation be harmonic 
and find 
@ = 


X |S 


(17) 


We shall show that the solution (15) constitutes the complete solution to our 
problem. To this end we change x; (i = 1, 2, 3) to x, y, z in (2) and using (17) 
write explicitly the components of stress acting on an element parallel to the 
plane z= 0: 

Tz, — (u,. a W, x) , 


Fy AO ie Wing) 3 (18) 
fs B 
2p (20,.—-— T). 


Since the partial derivatives of harmonic functions are themselves harmonic, 
and u, v, w and T are harmonic, this result shows that the stresses (18) are also 
harmonic. 

In order that the displacement (15) respresents the complete solution to 
our problem it has, first of all, to meet the boundary conditions on the plane 
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z =. These are in terms of stresses: 


,,.= 0,2, =0 -om 20. (19) 


Paes zy 


Furthermore, we shall require the fulfillment of the following conditions: 
a) the displacements and the stresses are regular except perhaps at the points 
of singularity, and except b) the components T,,, T,,, T,, Which are bounded 
everywhere in the region D + B; c) the resultant of the stresses acting on any 
surface enclosing a point of singularity, if any, and lying in D must vanish; 
d) the stresses must vanish in infinity”). 

Denote the stress vector acting on an element of area parallel to the plane 
z= const. and having the components T,,, 1,,, T,, by ¢,. It is bounded and 
vanishes on z = 0, by (19). Hence we may utilize Theorem 2 to conclude that it 
must vanish identically throughout the whole half-space. In other words, the 
stress field induced in the half-space by any steady-state temperature field is 
plane and parallel to the bounding plane. 

This theorem which follows as a natural consequence of the theory of har- 
monic functions extends the remarkable result of STERNBERG and MCDowELL [2] 
inferred in the special case of an arbitrarily prescribed surface temperature 
(bounded ‘region of exposure’). It is important to note in connection with the 
Equations (18) that although their right hand sides represent bounded functions, 
the boundedness does not concern the individual terms in (18) taken separately. 
Thus the derivatives of the components of displacement as well as the tem- 
perature may not be bounded. 

In the case of rotational symmetry around the z-axis we obtain instead of 
(18) by transformation to the cylindrical coordinates 7, z 


Tir =U (U,, 2 + wv, ‘i , | 


(18.1) 
t.,=(2w,,—£ Bi ; | 
where uw, is the radial component of the displacement. In this case we put 1,, 
w = $y, d3, the tangential component of the displacement vanishing identically. 
Again instead of (19) we obtain 
=i, = 0 on n=O. (20) 


Zr Zz 


Since the right-hand members of Equation (18) or (19) vanish identically 
throughout the region D + B, they establish important general relations be- 
tween the harmonic functions ¢; (¢ = 1, 2, 3). A trivial integration, neglecting 


18) The property d) is evident if one considers any hermispherical surface drawn in the body 
around a point on the boundary. For no resultant vector and moment of stresses acting on this 
surface, the stresses must go to zero with growing radius of the hemisphere. Apparently, we do not 
require that the displacements vanish at infinity and we shall see later that they may not vanish. 
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nonessential harmonic functions of integration, yields the following result in 
the non-symmetric case: 


u=b=— [ dsedz, 
v= $= 


w=$y— 3 [ Taz; 


| 

| 
— 
“ 
— 

= 
N 
cS 


(22) 


Denote 


®, = | Taz, (23) 


where ®, of course, is harmonic. A comparison of (21) and (22) with the cor- 
responding representation in [2]1*) shows that 


1 
eS ae (24) 


where y is the Boussinesq three-dimensional logarithmic potential applied in [2]. 

The strikingly simple result obtained in the preceding discussion permits us 
to solve a whole class of thermoelastic problems by a mere substitution of the 
expression for the temperature in (21) or (22), provided the temperature field 
for the prescribed steady-state heat flow is known. 

A detailed study of particular problems being beyond the scope of this 
paper, we now turn to a brief outline of some characteristic solutions connected 
with the half-space. 

1. At first consider a point source of heat located at the bounding plane 

z= 0. The solution of this problem was given by MELAN and Parkus [14], p. 74, 
using thermoelastic displacement potential and Love’s displacement function 
(the latter in order to remove the stress on the plane z = 0 created by the 
former function). The final result of [14] can be at once predicted from our 
previous argument, and the solution reduces to the following trivial reasoning. 

For a thermally insulated bounding plane 

De ae (25) 
where Q is the constant rate of heat introduced per unit time and K the heat 


16) Equations (31) and (32). 
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conductivity. By introduction of (25) into (22) we obtain the associated field 
of displacement 


ONT SP) oD 
bea  # wR REZ? (26) 
eee oe aE log(R + 2). | 


A simple inspection shows a singular behavior of displacements at the 
(singular) point of generation of heat, while at the point at infinity only w 
increases infinitely. The stress field generated by (26) renders the components 
t,, and t,, identically equal to zero, while the remaining nonvanishing compo- 
nents of stress become 


QOE«a il 
ez 4nKk R+2z’ | 
(27) 
QE« z 
Tod 4nK R(R+2)— 


These components tend to zero at infinity in accordance with the property 
d) required previously. 

2. The case of a given surface temperature in the region of exposure on 
z = 0 has been explored, as already mentioned, in [2] and [3], for two types of 
regions and distributions. Referring the reader to these papers for a detailed 
discussion, we borrow the formulae for the temperature field in the particular 
case in which a uniform temperature distribution exists over a circular region 
of exposure. 

In Sternberg-McDowell representation, we thus obtain 


T(r, 2) = 21K [Eth 8) — F(k, 0) — =] + BF, oh, (28) 
, 2 
E(k, 0) and F(k, 6) being the incomplete elliptic integrals of the first and second 
kind in the LEGENDRE’s normalform referred to the modulus # and the argument 
0, while K’ and E’ are the corresponding complete elliptic integrals referred 
to the complementary modulus k’. There designate 


Re oF k' = = (ar)? nm =[(a—7)? + 28", | (29) 
ro= [a + 7)? ame 
and 
i a? _v¥—a 
sin 9 = =e cos) = io (OPO S47) (30) 


a being the radius of the (circular) region of exposure. 


i 
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The representation by SNEDDON-LocKETT leads to the following alternate 
formula i in terms of the Bessel functions of the first kind: 


Fin.2)=T» f Te) Foleo) edo, (31) 


‘in which, in accordance with the notation of Eason, NOBLE and SNEDDON [15], 
_the improper integral can be also written down as J(1, 0; 0) with 


I (u,v; =f J hat) J,(0 t) ee Pdi (32) 


Since the elliptic integrals in (28) and the last integral have been extensively 
tabulated their numerical evaluation presents no difficulty. Of course, the two 
representations quoted above are equivalent to each other. This is shown in 
piS]. 
By utilizing either representation, e. g. (31), we obtain the associated field 
of displacement by a simple substitution into (22). In this way it may be 
trivially established that, in the case considered, 


u,=a(l+r)aTyI(1,1;-1), w=—a(l1+»)aT,I(1,0;—1), (33) 


in agreement with the result in [3], page 5. 

Since the determination of the displacement and stress fields in the half- 
space demands a mere knowledge of the associated temperature fields, many 
thermoelastic problems for the half-space may be at once solved provided the 
respective steady-state temperature problems are already solved. Let us quote 
here some more characteristic examples of known temperature fields which 
may be of interest: a) circular region of exposure at uniform temperature with 
no flow of heat over the rest of the bounding plane [16]; b) uniform flow of 
heat over a circular region with no flow of heat over the rest of the bounding 
plane [16]; c) uniform flow of heat over a circular region and zero temperature 
of the rest of the bounding plane [16]; d) a point source soaring over the boun- 
ding plane at a fixed distance [17]. 


5. The Infinite Thick Plate 


The foregoing analysis may be easily extended to the case of an elastic 
body bounded by two parallel planes, e.g. z= 0 and z=/h, but otherwise 
unlimited, i. e. to an infinite thick plate. Apparently, we might start with a 
theorem corresponding to the Theorem 21’). It is simpler, however, to choose a 
more heuristic approach. 


17) According to the oral communication of Professor C. Witcox, the proof of such a theorem 
for the infinite plate, sketched by him, involves subtle mathematical argument. 


142 Jerzy Nowi1nskt ZAMP 


Observe that since the stress vector f, on the planes parallel to the bounding 
plane of the half-space vanishes throughout the body, no interactions between 
the horizontal layers of the half-space exist. Thus we may separate a layer, 
say 0 <z SA, of the half-space from the remaining body without influencing 
the stresses in the disconnected parts. The complete solution to the problem is 
obtained by substitution into the Equations (21) or (22) of a solution of the 
equation of steady heat flow appropriate for an infinite plate. 

As already mentioned, the case of an infinite plate has been treated in 
some detail for a suitably chosen and axially symmetric surface temperature 
distribution e. g. in [2]. Hence we confine ourselves to a particular illustrative 
example of a point source of heat of strength Q located at the upper face of the 
plate while the lower face of the plate is kept at zero temperature. In this case 


Peo a Tae) sinh « (h — 2) 
ST | Sahu Pee ce 
0 


with K as thermal conductivity. 
By substituting (34) into Equations (22) we get in a trivial way 


Saf ie) Soke s 
f a : “COSHH Se 
0 
(35) 
,_ +49 f Jo(a7) cosh a (h — 2) 
oe 2K | a cosh « h ax, 


0 


in agreement with the result of LURIE [6], p. 196, obtained by a rather involved 


argument, using in the case of a layer equations similar to our general equations 
(21): 


6. The Sphere 


Let us now consider regions bounded by spherical surfaces. We shall confine 
ourselves to the solid sphere since the corresponding solutions for a thick 
spherical shell and a spherical cavity in an infinite medium may be obtained 
by a similar reasoning. : 

Take the center of the sphere at the origin of coordinates and denote the 


radius of the sphere by a. Let the sphere be immersed in an arbitrary steady- 
state temperature field!8). 


8 alte vam aS ; 

IN formal solution in integral representation for the sphere subjected to an arbitrary tem- 
perature field was obtained by BorcHarpr [18] in 1873. For the axisymmetric temperature field 
see [4] and Footnote 2. / 


‘ia 


wee eee ee 
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We assume the following expressions for the components of the displacement 


Nt ae) = 


4 eee mange | 
v=¢+ (PF —a)y,, ou, 
Be antes Ny | 


a 
; and we shall show that for a suitable choice of the functions ¢; (¢ = 1, 2, 3) 
_and »y they may represent a complete solution to the problem. 

We observe first that (36) becomes a solution of (5) if by virtue of (12) 


i ay ya iN [2 (1+ *) aT — 6) r4-1dr, (37) 
iy 
with the notation 
1—2 
00 = $12 + by + $32, A= Sees . (37.1) 


Hence the components of the associated field of stress become 
Ty, =2u (diet (*— a) y.,+2%y,+), 
Tyy = 2h (boy t+ (7? — 4) yyy +2yy,+®), 
t,2= 2 ($3,,+ (?—@)y.,+22y,,+9), 


Tey = [Ory + b22+2(7 — a) Pry t2V YY. + %Yy)] = 
Ty = [b12.+¢3,+2(—@) y,2+2(ey.+2Yy,)), 
Ty = lbo2+ b3,+2(7—@) py,+2(yy.t+zy,)], 
where 
DG, 21a oT. (39) 


In order that the stresses conform to the boundary conditions (4), in the 
absence of the surface tractions, the following equations must be satisfied by 
the functions ¢; (¢ = 1, 2, 3) and y on the surface 7 = a: 


2xd1, ¥ ba. + Zz Ps, x i v Pry “lp zy, . 
+2xry,+2ry,+0x=0, 

x by y = 2 be, a 2 os, +x ¢.,+ zo - 
+2yry,+2ry,+0y=0, 


XYi2t Vbn2t 2263. + 2930+ ¥ bay 
+2z2ry,+2ryp,+0z=0, 


(40) 
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with the notation 


d= 55 [Oy - (1 +9) a7). (41) 


1-2 


It may be shown by direct substitution that by virtue of Equations (5) the 


left-hand members of the above equations are harmonic. Thus since they 
vanish on the boundary = a of a bounded region they must vanish identically 
throughout the entire region occupied by the body, according to the Theorem 1. 
From this fact we derive the important conclusion that the only nonvanishing 
components of the stress vector ‘ (£, H, Z) which act on elements of spherical 
surfaces centered at the center of the body are 


te Fo HE 
= =2u Vv [Diego OV cg Mo ee 2) 
eee 
H=2p fo (Duy VY Pay X + Yyz z), (42) 
pene 
LZ = 2 — (0 222 + Pies ® Pad) - J 


We turn now to the system of equations (40) which we have found valid for 
the whole interval 0 <7 < a. Simple computation reduces the system to one 
equation for the function w 
IP 
Py, +2(L+yry,+ (L+%) pt+——-=0. (43) 


We assume now that the harmonic function p defined in the bounded simply 
connected region can be represented in a series of solid harmonics?*) 


v= S'a,(7)" Ys6, 0), (44) 


where a, are unknown coefficients and 6 and ¢ angles of latitude and longitude, 
respectively. We assume, furthermore, that the temperature, being an arbitrary 
harmonic function of coordinates, can be represented by a similar series 


T= S(2)" Y(0, 8) (45) 


Substituting of (44) and (45) into (43) and comparing of coefficients yields 
(1 + ») 


a, = — 
a 2[n? + (1+ 2) + (1+ »)] a (46) 
the series representation (44) being now complete. 
9) Evidently the solution of Equation (43) can be easily found in a closed form, but the method 


adopted in the sequel and modelled on the procedure of TREFFTz seems to be more convenient for 
practical purposes. 
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Our _next aim is to obtain the components t,,, t,, and 1, g of the stress 
vector #, resolved along the radius and the tangents to the meridian and the 
parallel at a point. By simple projection of (42) on the directions considered 


we find 
Tp, = bb (7? — a®) [2 (p,, cos?d + y,, sin? g) sin? 6 + 2p, cos? 0 
| + y.ysin2¢sin? 6+ (y,,cosd+y,.sind) sin 2 6], 


T¢ =- ie a a*) LP ex cos? d a Wiyy sin? d — Yizz 
+ yx, sin 2) sin26+ 2 (p,.cos¢d+y,,sin ¢) cos 2 6] , 


Tg = - (7? = a?) NG Y xx sin 2 d a DPiyy sin 2 ¢ 
+2y,,cos 2) sin 6 — 2 (y,,sind — p,, cos g) cos 8]. 


We observe that the computation of the components of stress acting on the 
elements of area perpendicular to the radii does not require the knowledge of the 
functions ¢; (¢ = 1, 2,3). This important conclusion greatly simplifies the 
following calculations. 
In order to find the effective equations for the stresses let us decompose (47) 
into the particular solutions corresponding to the separate terms of the series 
(44) and denoted by the superscript . In this way we find by lengthy but 
trivial computation : 


=I" (1 — S) (< i Y,, (9, 4) , 
wy =0(1- 5) (2) ¥,0 6.9). (48) 


=r (1-4) (<)" ao (8; )5 


a) sin 0 


where 


1 (1 +) (n—l)a 
ie (1+ 2+ (1+>») ~ ee) 


In the event that the arbitrary temperature distribution becomes sym- 
metric with regard to the axis through the poles, the surface harmonics reduce 
to zonal harmonics, or Legendre-polynomials P,,(p) with = cos 6. This reestab- 
lishes the known result for the axisymmetric problem derived in [4]. 

The preceding discussion hinges on the series representation of the function 
y (44). Since, by hypothesis, y is harmonic inside the region of the sphere its 
expansion in a series of solid harmonics is uniformly convergent in this region 
[13], p. 251 a. subs: Such a series, however, may also be regarded as a power 
series in three variables (x,y, z) according to the definition of the spherical 
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harmonics. Hence by virtue of the known theorem regarding power series”), 
any given partial deriivatives of y will be represented inside the sphere by a 
uniformly convergent series obtained by a termwise differentiation of the 
series (44). The series for the stresses T,,, Tg, Tg (48) being obtained in this 
way these series must be uniformly een acide the region considered. 
Since the aforementioned stresses vanish on the boundary 7 = a, we need not 
investigate the convergence there. 

Computation of the remaining components of the stress field Tyg, T,, and 
Tg as well as of the associated field of displacement requires the determination 
of the harmonic functions ¢; (i = 1, 2, 3) from (40). With no essential difficulty 
involved, this can be done following the procedure outlined by TREFFTz {11}, 
p. 104, for the isothermal case. 

It is obvious that using series representations (44) and (45) in negative or 
in negative as well as in positive powers of 7 we can obtain in a like manner the 
solution to the thermoelastic problems relative to arbitrary steady-state tem- 
perature distributions in a infinite medium exterior to a sphere (a spherical 
cavity in the infinite space) or in a thick spherical shell. 

_ In the foregoing discussion the components of the stress field have been 
represented in terms of surface harmonics which, bearing in mind the expansion 
(45), can be assumed as being known if the associated temperature field 
T(r, 6, 6) is known. Suppose for definiteness that the temperature distribution 
on the surface B of the sphere is prescribed, T(a, 0, 6) = S(6, ); S(6, d) being 
a sufficiently regular function of @ and ¢, say, its square is integrable over B. 


Then the series 
=)’ Y,(6, 6) (49) 


n=0 


converges in the mean to S(6, @) and the spherical harmonics Y,(0, 6) can be 
represented by the given function S(0, ¢) as well as by the extensively tabulated 
Legendre polynomials Pu) and the associated Legendre functions Pee 


where 4 = cos @. In fact, we have the following known representation: 


Y,(0; 0) = ol) H+ Ses )cos m @+ di”) s sin m ¢) P PM iu (50) 


m 


where for n = 1 


(B) 


Cn = Oma) | St S(0, 6) P\ (uw) cos m¢ do , A 


dM) = bm ‘i S(0, 6) P& (u) sin m ¢ do , 


(B) 


20) See e. g. [13], page 137. 
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and 

3 (2m+ 1) (n—m)! 

y Onin) aa »\ I (52) 
2a (n+m)!e,, 


Here e, a re Onn = 2if m=0 and P(u) = P,(u). For n = 0 we get 


Eimply Y, (8, 6) = c® acre c has to be computed from the first equation 
p (51). 


APPENDIX 


Proof of Theorem 2. Consider a harmonic function u(P) satisfying the 
_ requirements of the Theorem 2. 
‘ Let the half-space occupy the infinite region z > 0. To prove the theorem”), 
define u(P) for negative z as an odd function of z: 


u (x, y, — 2) = —u(x, y, z) . (a) 


First let us show that this new (extended) function u(x, y, z), — co <z< 0, 
is harmonic throughout the entire space. To this end, construct a sphere S 
with an arbitrary radius 9, and the boundary ’, the center of the sphere resting 
on the plane z = 0. 
Define a new function v(P) which is harmonic in S and which takes on 2 
_ the same values as u(P) on 2: 


v(P)=u(P) for Pond. (b) 


It may be easily shown that v(P) = 0 on the plane z= 0. In fact, an 
auxiliary function w(x, y, z) = 1/2 [v(x, y, z) + v(x, y, — 2)] is harmonic in S as 
a sum of two harmonic functions in S, and vanishes on »' because of (a) and (b). 
Hence, by virtue of the Theorem 1, certainly 


wP\=O for P on z=0. (c) 
But 
u(P) = w(P) for P on z=0 (d) 


by construction, hence the function v(P) vanishes on the plane z = 0. This 
plane divides the sphere S in two hemispheres. On the bounding surface of 
each of these hemispheres the values of v(P) coincide with the values of u(P). 
In fact, on the hemispherical surface this results from (b) and on the plane 
z = 0 both functions are zero, v(P) by (d) and u(P) by hypothesis. Hence, by 
virtue of the Theorem 1, the difference « — v vanishes inside any of the two 
eee eres and u = v throughout the sphere S. Hence the extended function 


21) The main line of the proof is modelled on SoBotov’s proof of the uniqueness of solution of 
the Dirichlet problem for the half-space [19]: In the argument we use the Liouville theorem which 
states that a harmonic function bounded in the entire space is constant. 
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u(P) is harmonic in S. Since the location of the origin and the radius of S are : 
arbitrary, this conclusion concerns also the entire space. 
Recall now that by hypothesis w(P) is bounded in the entire space. Hence 


by virtue of the Liouville Theorem, it is constant in the entire space. But it 


vanishes on z = 0, thus it may be only zero everywhere. This completes the 
proof of the Theorem 2. 
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4 Zusammenfassung 
Die Kelvin-Almansi-Darstellung einer Bipotentialfunktion vermittelst zweier 
7 harmonischer Funktionen wird zur Lésung elastischer Warmespannungsprobleme 
_ fiir stationare Temperaturfelder angewendet. 

| Es werden die Warmespannungen infolge willkiirlicher stationarer Temperatur- 
_ felder in einigen fundamentalen Bereichen, namlich im Halbraum, in der dicken 
_ Platte und in der Vollkugel untersucht. 


(Received: Juli 25, 1960.) 
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On the Dynamics of Turbulent Vortical Flow 


By ALan Reynocps, London, Great Britain?) 


The problem of the mathematical representation of the vortex flows occur- 
_ ring in nature has so far proved intractable. Two lines of attack have been used: 
the study of simple solutions of equations potentially capable of adequately 
representing the phenomena, and the study of flows with the necessary geo- 
metric complexity but satisfying simpler equations. The recent treatments of 
_ vortex flows by DEISSLER and PERLMUTTER, Lay, LoNG, PENGELLEY, RoTT, 
and VuLis and USTIMENKo ({1]?) to [6]) illustrate various approaches to the 
same problem. Every contribution represents a compromise between conflicting 
demands of reality; as the fluid properties are specified more realistically, the 
geometric verisimilitude must be correspondingly reduced. 

In most of the physical situations of interest the flows are turbulent. Never- 
theless, most investigators have found it expedient to base their studies of 
vortices in real fluids on the Navier-Stokes equations, extending their conclu- 
sions to turbulent flows by introducing a constant eddy viscosity to represent 
the activity of turbulence. 

Here the problem of the turbulent vortex will be attacked directly, using 
an order-of-magnitude analysis to pick out the important terms of REYNOLDs’ 
equations. This approach does not lead to explicit solutions unless further, more 
arbitrary assumptions are made. But it does simplify the problem and thus 
provides considerable insight into the mechanisms operating to set up the 
basic flow patterns. Also, the reduced equations would seem to provide a 
suitable basis for further attempts at the detailed prediction of vortex flows. 

From the outset we neglect the viscous terms in the momentum equations ; 
however, account will be taken of density variations, in order that the study 
be applicable to the interesting problem of the Ranque-Hilsch vortex tube. 


1) Imperial College of Science and Technology, Dpt. of Aeronautics. 
2) Numbers in brackets refer to References, p. 158. 
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} 
Thus the analysis can be valid only for flows in which there is a fairly high 
level of turbulent activity, high enough for the turbulent mixing stresses to 
far exceed the viscous terms. This neglect of viscosity means that the results” 
apply only in the turbulent body of the flow, away from solid boundaries where 
turbulence is suppressed. 

This investigation is particularly directed towards flows with large swirl. 
It is to such cases that the technique of approximation adopted is especially 
applicable. Thus the approximate results may be misleading if applied directly 
to ‘swirling pipe flow’ or other problems in which the swirl is small compared 
to the axial component of velocity. 

Finally, it should be noted that body forces are neglected. Consequently, 
only those flows in which the typical centripetal accelerations greatly exceed 
that of gravity are within the scope of this treatment. This stricture is likely 
to apply only to flows of liquids since .hese are commonly much slower than 
gas flows. 


1. The Governing Equations 


From continuity and EULER’s momentum equations we obtain, on taking 
mean values with respect to time and requiring the mean motion to be steady 
and axisymmetric: 


oO) — (4 a ee 
Pare )+ 5 (ew) = 90, | 
Ie Os7 6 v= O Fee 0 v2 oP 
y pine w) i 0 (9 uw) ey si sia Jay | 
1 
: —(rouwv) + iw) oe ee " 
: 2 gun + 2 @vw) + 2%" -0, | 
Ibe 8) = OF (esas oP 
mie (ro uw) + =e w?) = | 


whereo=R+ 0’, u=U+4', etc. We have then 


ouv=Ruv+o' (Uv'+Vw'4+u'v'), 
etc. 

The introduction of time-mean values has brought order into the system, 
but at a great price. For it is immediately apparent that this system of equations 
does not specify a determinate mathematical problem: we have now more 
dependent variables than governing equations. General principles on which 
further mathematical relationships might be based have not yet been devised. 
This is the fundamental problem facing us. 

As we cannot investigate turbulent vortex flows by obtaining and studying 
explicit solutions of the governing equations, we are forced to extract what 


4 
> 
- 
‘ 
, 
‘ 
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“information we can directly from the equations. To this end we shall simplify 
them by dropping the smaller terms, justifying these steps with an order-of- 
“magnitude analysis. Three cases will be studied: nearly plane motions, nearly 
helical motions as would occur within long tubes, and finally a more general 
class including these as special cases. 

2. Nearly Plane Swirling Flows 

Exact integrals can be found for some of the governing equations if axial 
variations are restricted. Indeed, these first integrals can be obtained even 
with the viscous terms retained if the coefficient of viscosity is assumed constant 
The results are applicable to flows that are the same in every section perpen- 
dicular to the axis. 

Consequently, in allowing small axial velocities and small axial gradients 
we are working away from a firm base. We can expect that neglect of less 
important terms will lead to considerable simplifications and that the accuracy 
of the approximate results will be reasonably high. 

The scheme for the order-of-magnitude analysis is as follows. We take 


0 0 
U,V~O), W~O@), ~ ~O(l), ~~ 06) 


to represent the basic features of the flow, and 
w?, vt, w®~ Of), gf ~ O(n) 
to represent the turbulence. 

Units are chosen such that 7, U, V, AP ~ O(1). Then R ~ O(1) also, for 
consistency of these assumptions with the equation governing tangential 
momentum; this result will be seen more clearly a posteriori. The quantities 
B, 6, €,n are, of course, assumed to be small compared with unity so that 
terms of higher orders in them can be neglected. 

Introducing these orders into equations (1) we find: 


rRU=C+ O(6e, BY? x) , 


oP 0 eg 
5 + yy RUY = —— + 068), 
2 
or eal ee 
oP iy Xa) — 1/2 4 
caae a a es R U w) — O(B 6, B ”) ’ | 
0(0) O(e, B) 


where only the lower order terms are indicated and those not marked ~ O(1). 
Note that residual errors due to non-zero compressibility remain in these 
relations even as 6,¢ > 0 to give purely planar flow. The axial momentum 
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relation appears to be the least accurate of the simplified equations; it has 
come into being as the result of introducing perturbations to the planar flow. 

Integration of the continuity and tangential relations has still been possible. 
A result of some interest can be derived from the integrated forms. We have 


PRu' ov =—G, for UV=0. 
If the radial variation of REYNOLDS’ stress does not change decisively with U, 


we can write, for small radial flows: 


& 


v2 


RUV=-— 


ioe as v= (se=") 1 


Cc i 

Thus it appears that V ~1/r for small or vanishing U's. This law is found to be 
obeyed in the fully turbulent region between two rotating cylinders, adjust- 
ments to particular boundary conditions taking place in boundary layers. But 
its justification for this case of zero radial flow is not simple unless the more 
general radial flow is considered first. 


3. Flow in a Vortex Tube 


Here the system of approximation is 


V~O), U~0~), 2 ~o), 2 ~ 06), 


u'2, v’, w’2 ~ O(B), 0’ ~ O(n). 


Units are again chosen so that r, V, 4P ~ O(1); as before, R ~ O(1), for 
consistency. 

The geometric specializations imposed seem natural in a problem of high- 
swirl flow in a tube, although they will not necessarily be valid near the ends 
of the tube if there are constrictions there. Introducing them into the continuity 
equation (1), we find that 


AT p é 
W~ Ol, ), 
at most, so that 
Ie. 0 , ‘ 
5 RU) +2. (RW) = OBI? n). 
O(e) O(e) 
Also, 
oP ae 
va == R y + O() ’ 
UR VW) + see gall (7? Ruv) = O(8 6, Be 
Oz Ole) v2 Ov Bes tae : B ") : (3) 


OP ee ae ed #3 
oz + ag (RW?) + —- 5 Raw) = 0(8.8, | fin). 
0(6) O(e*/6) O(B, €*/6) 


a ate 
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The accuracy of these simplified results does not seem to be as high as that 


of equations (2). However, detailed investigations of the flow within a vortex 
tube have been made by the present writer; these experimental data allow 
estimates of the errors to be made for this particular case. 


The vortex tube was about 7-5 cm (3 inches) in diameter and about 120 cm 


| (48 inches) long. Air entered through slots cut tangentially through the wall 


near one end and escaped to the atmosphere through throttling devices at the 
ends of the tube. Internal flow measurements were made with a variety of 


blockages at the ends of the tube — solid end plates, orifice plates, perforated 
plates, or no blockage at all. Eight configurations were studied for overall 


pressure ratios of 2-0 and 2:3. 
In the turbulent body of the flow we find, typically, W/V ~ 1/4, suggesting 
that « ~ 6/4. Considering the pressure distributions we find, typically, 


AP A 
Ar 50m 


relating changes down the length and across the radius of the tube, suggesting 


that 6 ~ 1/50, e ~ 1/200. The experimental results indicate that for balance 


of the tangential and axial momentum equation we must have f not very dif- 
ferent from 6 and e. Say B ~2¢ ~ 6/2 ~ 1/100. This implies that the intensity 
of turbulence in the tube ~ 10%. As this is of the order erecieg, it offers a 


slight check on these estimates. 

Finally we must estimate 7, the order of the turbulent density fluctuations. 
The following considerations allow us to do this. 

If the fluid is homogeneous (that is, such that it would be homogeneous on 
being brought to rest) the density fluctuations will be dependent on the turbulent 
velocity intensity. It seems plausible that these density differences are main- 
tained by the centripetal pressure gradients within the eddies forming the 

turbulent velocity field. The pressure difference from core to periphery of 
a ‘typical’ eddy is related to the velocity field by 


p' ~ Ru? ~ Off). 


We know that 0’/p’~ K, the compressibility of the fluid. Hence, 4 ~ K B. 
We have now only to estimate K, the compressibility, in the units appro- 
priate to this study. We have taken AP, the pressure difference across the 
vortex, ~ O(1). But in the high velocity air flow of the experimental situation 
we have also AR, the density difference across the vortex, ~ R ~ O(1). Then 


AR 


“K~ 55 ~ 0(1) 


in the appropriate units. Binally, ny ~ 1/100. 
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We can now give estimates of the error terms in this situation as percen- 
tages of those retained: 


: lt RU) + (ek W) = 0(20%) 1 term | 
ie mais — + O(1%) , 2 terms 
ats V W) ae a oe R u v) a O(15%) , 1 term t (3a) 
oF + (RWWA tS Raw) 02%), “lien 
or 
a = el Ru’ w') = (10%) , 2 terms. ) 


The number of the larger rejected terms is indicated to the right. 
Very likely the errors estimated for the continuity_and tangential momen- 
tum equations are much too high. The more accurate forms of these equations 


are 


Tee 0 Log Pita som oe 
~s- RU) + (RW) + — —(ro'u') = 0(6 Bi? 0), 


iL) 


1 
v2 cal 


P Rud) + sees VW) + 1, 2 (PV 0'w) = Oe, Bd). 


y2 


The estimates of the mixing terms were made.on the basis of 100% corre- 


lation; for example, ou! ~ O(61? y). In reality, the correlations will be much 
lower, and the error terms dependent on them correspondingly reduced. It 


seems likely that 0’ wu’ > 0 in vortex flows, since lighter lumps of the turbulent 
fluid will tend inwards under the action of Archimedean buoyancy forces. This 
finite correlation will definitely introduce errors into the simplified equations, 
although much smaller errors than those estimated above. 

Finally, we conclude that the terms neglected in the simplified equations 
are probably only a few per cent of those retained. The simplified forms should 
represent the mass and momentum balances with accuracy sufficient for most 
purposes. 

Although it is not possible to obtain integrals for this case, some of the 
gross features of the flows in the central part of the vortex tube can be grasped 
using the approximate relations (3, 3a). 

The tangential momentum relation is no longer a simple balancing of the 
Reynolds’ stress term against convection of mean vorticity; the effects o! 
axial convection are not negligible here: 

RU (|e te =| ate Rw oe ys re (P Ru'v'). 


y* Or 
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(Note that the simple continuity relations are sufficiently accurate in both 


2 systems of approximation to allow the mean-motion acceleration terms to be 


_ written in their more usual forms without altering the orders of the errors.) 


7 


The simple radial equilibrium relation, 


retains its applicability in this type of flow. 


We have found in equations (3a) that the axial momentum equation can 


_ be simplified even farther, by dropping the mean-motion accelerations, without 
- introducing disastrous errors. Although the final result is the least accurate 
of the set, it us sufficient to show that axial pressure changes are chiefly 
- accountable to REYNOLDs’ stresses. This result is in accord with the conclusion 
_ reached upon examination of the experimental results on which the estimates 
_ of the errors are based. 


4. The Vortex Shape Hypothesis 


A key feature of the two systems of approximation just considered is the 
relationship between the gradients within the vortex, the requirement that 


- changes of any mean-value function be much more gradual in the axial direction 
_ than in the radial direction. It seems appropriate to term this a vortex shape 
hypothesis. The considerable simplifications of the equations of motion are 
chiefly accountable to the introduction of this concept. 


The hypothesis has needed little justification in the two cases that have been 
examined. In the former, small perturbations were imposed on planar flows; 


in the latter, the elongated chamber in which the motion took place offered 
_ good reason for the assumption. But, as we shall now see, this approximation is 


useful for other vortex problems to which its applicability is not so immediately 


apparent. 


On examining the swirling flow of water through a contracting conical 
chamber, BINNIE and TEARE [7] found a cylindrical air core. This parallel- 
sided streamsurface indicates that the hypothesis is valid away from the 
tapering sides of the chamber. The study of swirling flow through a convergent- 
divergent nozzle (BINNIE, Hookines, and KAMEL [8]) revealed again a cylin- 
drical air core in the contracting subcritical flow above the throat. The super- 
critical flow below the throat did not form a cylindrical core along the axis but 
fanned along the diverging walls. However, a supercritical flow will not neces- 


sarily violate the assumption that radial changes are very much more rapid 


than axial. If the flow is retained within a cylinder, the hypothesis may still 
be valid. 


. 
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A somewhat similar situation is found in the vortices shed from lifting 
surfaces. Commonly they tighten in on the vortex core and the elongation 
hypothesis represents the pattern well. But a vortex burst phenomenon is 
sometimes encountered in trailing vortices [9], the tight core breaking up and 
the swirling fluid moving away from the axis. Another case of interest is a 
swirling jet emerging from an orifice or pipe into the atmosphere. Commonly” 
the vortex spreads out, often leaving the core free for reversed flow. 

Finally : there are cases for which the vortex shape hypothesis is not useful, 
but is does adequately represent the situation in the heart of the vortex for a 
variety of problems, not only flows confined within nozzles and tubes, but also 
some cases of unrestrained flow in the atmosphere. 

The reasons for its wide applicability and the limitations to its applicability 
are not hard to see. In the first place, the radial and swirl velocity components 
must be zero on the vortex axis. Therefore, their gradients must be related as 
suggested by the elongation hypothesis. If, further, we have rapidly swirling 
fluid close to the axis (that is, a high angular velocity) the radial pressure 
gradient will be high in the core. This in turn leads to steep radial density 
gradients, whether by air core formation in a liquid or by the normal effect of 
compressibility in a gas. 

Thus we see that the kinematical requirement of axisymmetry, together 
with the dynamical requirement of radial equilibrium for rapidly rotating flow, 
must specify a flow in which radial gradients are large compared to axial 
gradients. It is to be emphasized that this relationship will hold for the gra- 
dients of all the flow variables only if there is a high angular velocity in the fluid. 

It is interesting to note that many of the highly developed branches of fluid . 
mechanics depend on the assumption that changes in one direction are small 
compared to those in another. As examples we can mention boundary layer 
theory, the classical theory of surface waves, and the theory of the thin wing. 
The power of this hypothesis is not surprising. For fluid mechanical theory 
thrives on simplification of its initially difficult equations. And in a problem 
governed by equations with two independent variables a great simplification 
can almost always be made if the gradients in one direction are small. It seems 
safe to predict that this hypothesis has not yet reached the limit of its use- 
fulness. 


5. General Flow with High Swirl 


We have seen that the vortex shape hypothesis used previously for the two 
special cases has a wider validity in high-swirl flows. We note also that the 
two sets of simplified equations (2, 3) contain, for the most part, the same 
terms. That is, only a few terms become important or lose their importance as 
we shift from one sort of flow to the other, while most of the terms of REy- 
NOLDs’ equations (1) are negligible for both categories. 
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It seems plausible that by combining the two sets, retaining the dominant 
terms of each class, we will obtain equations valid for a less restricted group of 
vortical flows, a group containing not only the two specified classes, but also 
the whole range of flows intermediate to them. 

The combined set is 


Sy RU) + 


or aah or 


(RU), 
(4) 
(Ruy) =0, 


7) 

0z 
0 
or 


OP ait 8) ne 
= oy" Ruw) =0. 


We know that these equations are valid for large swirl so long as one of the 
other velocity components is small. But are they also valid when all the velocity 
components are of the same order? This is easily checked. Taking 


“ ~O(6), wu’, v'2, m'? ~ O(8), 


U,V.W~ol), 2% ~0(1), 


we obtain the approximate forms 


+RU=C 2 Ol, 


We v 1 2 
op = RS Wy wT RU) + 06,8), iS) 


PRUV=G+ 068), 
ryRUW =F + 0(6, 8). 


These results are contained in the general set (4). For this degenerate case 
all the REYNOLDS’ stresses are found to be negligible and explicit solutions are 
possible. 

The problem of the turbulent vortex is still intractable in the simplified 
form (4): we have two more dependent variables than governing equations. 
However, the order-of-magnitude analysis has been useful in displaying the 
essentials of the mechanics of flows with high swirl. 

The effects of compressibility are now represented solely by the time-mean 
density; the continuity and radial equilibrium relations are not directly in- 
fluenced by turbulence; only two of Reynolds’ apparent stresses are found to 
be important in determining the motion — those acting on cylindrical surfaces. 
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Notation 
Y radial co-ordinate; 
Zz axial co-ordinate; 
u, U,w’ radial velocity component, its time-mean value, and the turbu- 
lent perturbation to it; 
OVa oe tangential or swirl velocity component, etc. ; 
w, W, w’ axial velocity component, etc. ; 
DB, Pye. fluid pressure, etc. ; 
0, R, 0’ fluid density, etc. ; 
6, 6, €, 7 numbers small with respect to unity; 
C,F,G,G, constants of integration. 
REFERENCES 
[1] R. G. DeissLter and M. PeRLMuTTER, An Analysis of the Energy Separation in 


[2] 


[9] 


Laminar and Turbulent Compressible Vortex Flows, Proc. 1958. Heat Transfer 
and Fluid Mechanics Institute, Stanford University, Palo Alto, California. 

J. E. Lay, An Experimental and Analytical Study of Vortex-Flow Temperature 
Separation by Superposition of Spiral and Axial Flows, Trans. ASME, 87 [C], 
3), AWA, Zils) (UL) 

R. R. Lone, Sources and Sinks at the Axis of a Rotating Liquid, Johns Hopkins 
University, Civil Engineering Dept., Techn. Rep. 9 (1957). 

C. D. PENGELLEY, Flow in a Viscous Vortex, J. Appl. Phys., 28, 1, 86 (1957). 


] N. Rott, On the Viscous Core of a Line Vortex, Z. angew. Math. Phys. IXb, 5/6, 


543 (1958); X, 1, 73 (1959). 


| L. A. Vuts and B. P. Ustrmenko, Uber die Aerodynamik der Zyklonfeuerungs- 


kammern, Energietechnik, 5, 6, 265 (1955). 

A.M. Binnie and J. D. TEARE, Experiments on the Flow of Swirling Water 
through a Pressure Nozzle and an Open Trumpet, Proc. Roy. Soc. Lond., A 235, 
78 (1956). 

A.M. Binnigz, G. A. Hooxines, and M. Y.M. Kamer, The Flow of Swirling 
Water through a Convergent-Divergent Nozzle, J. Fl. Mech., 3, 261 (1957-8). 

B. J. ELLE, On the Breakdown at High Incidences of the Leading Edge Vortices on 
Delta Wings, J. Roy. Aero. Soc., 64, 596, 491 (1960). 


Zusammenfassung 


Die Wirbelstr6mung in einer turbulenten kompressiblen Fliissigkeit wird behan- 


delt. Die Bewegungsgleichungen werden durch die Vernachlassigung der kleineren 
Glieder vereinfacht. Die experimentellen Ergebnisse erlauben eine Fehlerrechnung 
fur einen Fall. Ein wichtiger Schritt in der Vereinfachung ist die Voraussetzung, 
dass die Radialgradienten viel grésser sind als die Axialgradienten. Fiir die kraf- 
tigen Wirbel erwies sich diese Annahme sowohl durch die Erfahrung als auch aus 
den wesentlichen theoretischen Betrachtungen heraus als begriindet. Die verein- 
fachten Gleichungen zeigen einige wichtige Aspekte der schnelldrehenden Stré- 
mungen und geben einen Anhaltspunkt fiir die ausfiihrlicheren Studien. 
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An Investigation into the Mechanism of Ice Crystal Nucleation 
by Proton Spin Resonance Spectroscopy 


By GEOFFREY T. BARNEs and RaymMuND SANGER, Ziirich!) 


In experiments on weather modification silver iodide is dispersed in the atmo- 
sphere and under certain conditions causes marked changes in the pattern of 
precipitation. These phenomena are due to the activity of silver iodide as an ice 
forming agent, but the mechanism of this process is not yet fully understood. This 
paper reports on some preliminary experiments in a program aimed at elucidating 
the mechanism of ice nucleation by silver iodide and other nucleants. 

Three mechanisms have been considered in order to account for the formation 
of ice crystals when nuclei are introduced into a cloud of supercooled water droplets 
(e. g. see Mason and VAN DEN HEvVEL [1]?) and SANGER [2]). 

(a) In ‘freezing’ nucleation the nuclei are thought to act by collision with 
supercooled droplets in the cloud, causing them to freeze. Statistical considerations 
and several experimental observations (Katz [3] and KoeEnic [4]) indicate that this 
process is relatively unimportant in cloud seeding experiments. 

(b) In the ‘pure sublimation’ process the nuclei adsorb water vapour from the 
atmosphere directly into an ice-like structure. 

(c) In what might be called condensation-freezing nucleation water vapour is 
first adsorbed onto the nuclei to form a liquid layer which later freezes. 

In all three mechanisms the embryo ice crystals thus formed can grow further 
by transport (through the vapour) of water molecules from nearby liquid droplets 
in the cloud. 

All three mechanisms are idealized cases as they assume that the surfaces of 
the nuclei are initially free of adsorbed water, but, unless special efforts have been 
made to achieve this condition, the assumption is undoubtedly incorrect for most, 
if not all, nucleants. For example, BiRSTEIN [5] has shown that at temperatures 
as high as 20°C and at 50% relative humidity there are some four to eight molecular 
layers of water on AgI and PbI, powders. We must therefore assume that, normally, 
ice forming nucleation follows a mechanism in which an appreciable amount of the 
water required to initiate the growth of ice crystals is already present on the nuclei 
when they are injected into the cloud. 

The activity of the introduced nuclei thus depends on the physical state of this 
adsorbed water when the particles have assumed the temperature of the cloud. If 
the adsorbed water is all in a liquid-like or disordered state there will be no activity ; 
if some or all is in an ice-like or ordered state the nuclei may be active if other 
factors, such as the aqueous vapour pressure, are favourable. Since particular 
nuclei may be inactive at one temperature and active at some lower temperature, 
these ideas imply that there must be a phase transition involving at least some 
of the adsorbed water as the nuclei change from an inactive to an active condition 
(or vice versa). Therefore, in order to demonstrate the feasibility of this ice nucleation 


1) Laboratory of Atmospheric Physics, Swiss Federal Institute of Technology. 
2) Numbers in brackets refer to references, page 163. 
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process it is first necessary to show that such phase transitions do occur on ice 
forming nuclei. This is the purpose of the experiments described here. In some 
other systems transitions in adsorbed layers have been observed previously (e. g. 
FIsHER and Mc MItran [6]; Mays and Brapy [7]. ; 

It should be noted that the adsorbed water need not (and usually will not) be 
uniformly distributed over the surface. Nor is it necessary for the transition of the 
whole adsorbed layer on a particle to occur at the one temperature: water films 
on different parts of the particle may well have quite different transition points, 
either because of variations in the thickness of the adsorbed film or because of 
differences in the surface structure of the nucleus. Furthermore, the proposed model 
does not rule out the possibility that the first adsorbed layer is more or less im- 
mobile even at temperatures where the overlying layers are fluid. 

At first sight it would appear that the formation of an ice embryo on a parti- 
cular nucleus should be a function of the temperature only. However in experiments 
with silver iodide in clouds formed from salt solutions of various concentrations 
Katz [8] has observed that a particular activity value*) was always reached at the 


same degree of supersaturation with respect to ice and not at the same temperature. _ 


Two alternative explanations consistent with our hypothesis on ice forming 
nucleation can be advanced: firstly, that the degree of supersaturation has no 
effect on the formation of stable ice embryos, but does markedly affect their 
further growth into visible ice crystals; or secondly, that initially the ordered 
zones are unstable and will revert to a disordered condition unless the supply of 
water vapour is sufficient to cause rapid growth to a stable size. 

The model suggested would also provide an explanation for the effect of particle 
size, as observed by HosLEeR and SPALDING [9] and by SANO, FuJITANT, and MAENA 
[10]. 

To demonstrate the existence of phase transitions in the adsorbed water films 
on ice forming nuclei we have made use of the nuclear spin resonance of the protons 
in the adsorbed water molecules. If the molecules are able to move fairly freely the 
local magnetic fields due to neighbouring water molecules tend to average out and 
the nuclear magnetic resonance (NMR) spectrum exhibits a sharp proton resonance 
line. However if movement is appreciably hindered, as in an ice-like structure, this 
averaging out is no longer possible and the local fields will vary somewhat from 
proton to proton leading to a broadening of the resonance line. 

The NMR spectrometer is basically that described by Primas and GUNTHARD 
[11] modified by Arnpr [12] for low resolution work with solids. The B, field is 
provided by a permanent magnet of approximately 6,000 gauss, corresponding to 
a proton resonance frequency of 25 megacycles. The usual modulation method for 
solid state NMR was used, giving the first derivative of the absorption curve. In 
the present experiments the modulation amplitude was 43-2 milligauss. With AgI 
the sweep rate was 2:41 _m gauss/sec, and with CuO it was 4:88 m gauss/sec. In 
both cases the time constant of the final amplifier stage was set at 2 sec. 

Samples of AgI, CuO, Cu,O, and CuS from the same batches as those used by 
Katz [8, 13] were prepared by allowing them to stand in air saturated with water 
vapour at 0°C for several days. The surface areas of these samples have been 
determined by ZETTLEMOYER [14]; the values found by water vapour adsorption 
are given in Table 1. 

Examination of the samples at room temperature showed a well defined proton 
resonance signal with CuO, a weaker signal with AgI, and no observable signals 


8) The activity is defined as the ratio of the number of ice crystals observed in a supercooled 
cloud to the total number of nucleating particles introduced into the cloud. 
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with Cu,O or CuS. Accordingly only AgI and CuO could be used in studying the 
effects of temperature on the proton resonance signal. Owing to the low concen- 
trations of protons in these samples (ca. 10!® protons within the detector coil 
‘compared with 10 for liquid water) it was necessary to operate the spectrometer 
near its maximum sensitivity. Consequently the signal to noise ratio was rather 
low and precise delineation of the spectra was not possible. Figure 1 shows some 
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Figure | 


Typical proton spin resonance spectra (including noise) of the adsorbed water films on Agl and 
CuO powders at various temperatures. 


typical records. Obviously with records similar to that of AgI at — 15-5°C (Figure 1) 
the existence of a spectral line can only be ascertained after scanning the region 
several times. 

When the sample temperature was lowered there was at first no effect on the 
resonance signal, then at a certain temperature the wave height began to decrease. 
However the line width remained approximately constant throughout the entire 
temperature range with both samples, indicating that with decreasing temperature 
the observed signal from liquid-like adsorbed water was being replaced by a much 
broader, shallower signal from the ice-like state which we were unable to detect. 

The effects of temperature on the amplitudes of the resonance signals (in arbi- 
trary units) are shown in Figure 2. The wave height data for the two different 
samples cannot be directly compared as different B, (exciting) fields were used 
and no precise calibration was available. However this is of little importance at 
present as the most interesting observation for our purpose is in the comparison 
along the temperature scale. 

For both samples the curves exhibit breaks which may be interpreted as the 
ynset of ordering or ‘freezing’ (with falling temperature) or the completion of 


ZAMP XII/11 


162 Kurze Mitteilungen — Brief Reports — Communications bréves ZAMP. 
disordering or ‘melting’ (with rising temperature) in the adsorbed water films. It 
is interesting to note also that the change in wave height, and hence the order- 
disorder transition, was not abrupt but was spread over a wide temperature range 
of more than twenty degrees. This spread could be due to differences in the 
transition temperatures for different particles in the sample or to heterogeneity 
in the films on individual particles. Probably both factors contribute. 
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Figure 2 
The effect of sample temperature on the wave height of the proton spin resonance signal from 
water adsorbed on AgI und CuO powders. 


Thus it is clear from Figure 2 that the adsorbed film on AgI has a much higher 
transition temperature than that on CuO, and in accordance with the model for ice 
crystal nucleation outlined above there should be a corresponding difference in the 
temperatures at which these powders act as nuclei. This relationship is shown in 
Table 1, using the nucleation data of Karz [8, 13]. 

In comparing the results in Table 1 a number of complicating factors must be 
considered. Firstly, the high noise level in the NMR spectra introduces considerable 
error in the measurement of wave height; this is especially noticeable in the curve 
for AgI with decreasing temperature. Secondly, whereas it is possible to detect one 
active nucleus in 10% in the cloud chamber experiments, it is impossible to observe 
a decrease in wave height with this sensitivity. Thirdly, other conditions, such as 
the aqueous vapour pressure, may influence one or both types of experiment. 

Bearing these difficulties in mind, the agreement between the nucleation and 
NMR data does appear to be good enough to support the mechanism of ice forming 
nucleation which we have discussed above. Further work with substances of higher 
specific area should yield more accurate spin resonance data and enable a more 
quantitative appraisal of the suggested nucleation mechanism to be made. 
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Table 1 


Specific surface areas [14], ice crystal nucleation data [8, 13], and the onset of ordering or 
freezing’ (transition temperature) in the adsorbed water films on nuclei as determined 
; by proton spin resonance measurements. 
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Surface Area nea Nucleation Transition 
Substance (by HO ads.) (fae pie he Temperature Temperature 
if ° ° 
(m*/e) active nuclei) ee) ee) 
ee et et a) 
ee —4 +5(+5) (cooling) 
Agl 0-1 g 
8 . 10-3 ==6 —1(+3) (warming) 
104 —18 : 
CuO 4-7 ae tees —22(+3) (cooling) | 
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Zusammenfassung 


Unter Heranziehung der Kernresonanzspektroskopie wird versucht, den Vor- 
yang der Keimbildung an Eiskeimbildungskernen wie AgI und CuO abzuklaren. Die 
Ergebnisse der ersten Versuche zeigen, dass die adsorbierte Wasserschicht unter- 
1alb einer bestimmten Temperatur, die nahe der Schwellentemperatur fiir die Eis- 
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keimbildungsaktivitat liegt, teilweise in einem fliissigkeitsartigen und teilweise in 
festkérperartigem Zustand auftritt. Bei abnehmender Temperatur wachst der feste 
Anteil der Schicht auf Kosten des fliissigen Bestandes. Der volle Verfestigungs-_ 
vorgang wickelt sich dabei in einem Temperaturbereich von ca. 20°C ab. Die 
Beobachtungen legen es nahe, anzunehmen, dass zwischen dem Auftreten der 
eisartigen Struktur in der adsorbierten Schicht und der Eisbildungswirksamkeit — 
der Kernstoffe eine Beziehung besteht. 


(Received: January 3, 1961.) 


Two-Dimensional Steady Temperature Fields 
in a Stratiform Half-Space 


By VActav VopiéKa, Plzen, Czechoslovakia 


1. Statement of the problem. The half-space 
42K, —Oaya to, = cw—%7 = + co 
is made of m homogeneous isotropic parts 
Mpeg ss MSS Wp, * = CO <= Y= co, = cones = co, 1S hk See 
WW a OOS YS a OS, 0 = ss = 00 


with thermal conductivities 4,, respectively. The question is to determine the 
steady temperature distribution in the solid, provided that the front face x = x4 
is kept at the temperature /(y) and that there is continuity both of the temperature 
and of the flux at the surfaces of separation between several parts of the body. 

Mathematically speaking we have to find the solution uw, = u,(x, y) of the 
equations 


O°u, | OPuy 


Ox?! Oy® a Ae ay SM SING OY Se ay OO, Salis hE t caln eely 
07u,, O7U, : : 
Ox? Oy = 0, 4 > Xy_p ee ah Se A See, 2 (1.2) 
with the following conditions: 
U(%, ¥) =f(v), -co<y<+o, 
(2.1) 
Un(Vp V) = Unyil% eV), — 00 <y< + oo, 
OUp OU yp 4. 
Ay = — itl x=, —oocye toc, A, ae <ksn-1 
Ox OA f ; ’ x Us terri gt ay Le 
(202) 


jim n(x, y) = 0, — 0 <y< +400. (2. 3) 


x 
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2. The Fourier transforms of u,(¥, y). Making the usual assumptions on the 


disappearance of wu, and oe for |y| > oo the familiar transformation 
- at oo 
= 1 pee 
Hn) = FU) = —— | my)” ay, 
V22 
—oo 
(3) 
ee 
ny) = F* (h(n)] = —— | Tin) ean | 
y2 Fp 


changes our problem (1. 1)—(2. 3) into the equations 


Pu, 


xt ae Up = Oe Nyy a Vip k= ai Ze oo gp (RS Ale, (4.1) 
au By 

ORE Ahan 0 ® Heh (4.2) 
U1(%, 9) = f(n) (5.1) 

y Z di, di, 
Up 41 (Xp 1) = U,(4p, ») , — Ay F ae = 44> Vhs Rk => i, a Ano 4d diss (5. 2) 
i U, = SNS) 
Poe Un(*, 17) (5. 3) 


for determining the Fourier transforms u;,(¥%, ) of the required functions u,(4, Vv). 
Integrating (4.1) and (4.2) gives generally 


U,,(%, 9) = A; cosh (* — *,_,) 7+ B, sinh (¥—4,_,)7, 1 Sk Sn, (6) 
where the coefficients A,, B;, os to be found from (5.1)—(5. 3) i. e. from the equations 
A; =f(n); A,+ 2B, =0, (7.1) 
Ay i, = A, coshl, 7 + B, sinh hl, 7, 
Buys = Ay (Ay Sinh Gn + B, cosh tl, ), ly = %p—-%e_-1 1 Sk Sn—1. (7.2) 
Introducing the matrices 


lp, 
B,, | 


| cosh/,7, sinh J, | ‘ 
| , 


Sksn-1, 


Ky = > Ll Sk S72, My(n) =| 


A, sinh J, 7, A, cosh 1, 4 | 
Qin) = M;(n) My_a(n) -** Ma(n) Man), 1 Sk Swi, (8) 


the well-known manipulation with (7.2) and with the first equation (7.1) gives 


us [2]*) = Sot 
Kaya = Onin) |? l ee ont (9) 


I 1 |i 


1) Numbers in brackets refer to References, page 168. 
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(1, 9 ( (7) 


, 


Writing Q,_ = = pee y,s =1, 2, expressing A, and B, from the last 


equation (9) oral ee into the second condition (7.1) yields the value of Beg 
On the whole, one finds the following expressions for the coefficients A,, B,: 


A= fn), Br= — yey tea), 
Axirl] __ Qeln) | NO) la, i 
ie N(x) 2h | LSA an —" i, 


Z(n) = a) +O), No) = 0) + oS, (10) 


7(n) is the Fourier transform of f(y). 
The main difficulty lies in finding the matrices Q,(7). The first of them, Q,(y), 


Q.(7) and Q3(n), are.given explicitely in [2]. 

3. General solution of the problem (1.1)—(2.3). Finding A; = A,(m), B, = B;(n) 
from (10) and substituting into (6) gives u,(¥, 7) and from this one obtains, by aid 
of (3), the required solution 


+co 
U;,(*, Y) = cant [A, cosh (¥ — ¥,_,) 7 + B, sinh (w — %,_,) nJe —iny dn, 
Lleken See 
4. A homogeneous half-space. In this case we have A, = A, =---=A,_,=1 


and the matrices Q,(y) are expressible in closed form 


|| cosh (4, — %)) 7, sinh (x, — x,) 1 
x(n) = | k o) 7] (¥4 0) 77 | 


5 USES Sl. 


sinh (4, — %9) 7, cosh (%, — x» )n 
Using this, the formulae as give 
(n) = N(n) = e(*n-a~*y) ", 
A,=— B,= jin), 
Aig = = Beg ee RT Oi) Lae = wed 


and the solution (11) appears in the form 
U2; 4) == ———— / f(n) e~ © —% +)" dy (12) 
TU « 


which is valid for the whole extent of the solid. 
Using the well-known result 


ere Tee bie 

Soin ae COS 7 @ eR Tl OO 

4 es i] a? + 02 dg = = / ‘a2 + 02 dg 
0 @ 


of Laplace shows ~ see (3) — that 


2 1 
=O Eee ae eee. 
e -«/2 laren (13) 


= 
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_and the convolution theorem for Fourier transforms changes the solution (12) into 
its customary form [3] 


+00 
Add! Hone ve 
ee ela i (v= x)? + — a)? 7 


5. A half-space of two parts. Here we have » = 2 and the relations (10) give 
Z(n) = e*"+4 (Ay—1)sinhi,y, N() = e+ (A, — 1) cosh], 7, 
_N(n) [A, cosh (x — x9) 1 + B, sinh (x — x9) 7] 
= f(n) [e~ )" + (Ay — 1) cosh (¥, — x) ni), 
N(n) [A, cosh (¥ — %,) n + B, sinh (¥ — x,) 9] = Aye — 7 7 f(n) 


_and the solution (11) appears in the form 


; pect 
1 + f(n) (4,4) — tny 
Uy(¥, y) = ——— fF -~' [e™ 74 (A, — 1) cosh (4, — x) yn] e "dn, 
1 y) y2a- N(n) ( Z ) ( a ) ] 1] 
cre 
= Ay f() (x, — x) n — iny 
Ux(%, ¥) = ix / Nin) é dn , 
N(y) = e(~ 7) 4 (A, — 1) cosh (*, — %)-7. (15) 


If especially A, = 1, then our result (15) changes into (14). Generally, the inte- 
grations in (15) are to be accomplished by numerical methods, for the exact cal- 
culations are exceedingly difficult. Apart from the above special case A, = 1 and 
from the artificial cases characterized by A, = 0 or A, = oo the author of the present 
paper does not know any situation with arbitrary f(y), where it would be possible 
to express the solution in any kind of closed form like that of ref. [2]. 

Nevertheless, there exist special functions f(y) admitting fargoing simplifications 
of the above complicated result (15). To give an example let us take 


1,2 - 
fy) =-—"" : c< y<+oo. (16) 
Or Yah +) 


Writing this in the form, 


and using (13) gives 


am (A, — 1)/'a ein j (Ay +1) x ein 
Yorn rez: V2 


= — f(A, — 1) e- "74 (A, + 1) eo") = 2274 (A, — 1) coshd, 9 = N(x). 
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Hence, applying again (13), the solution (15) becomes 


we WE, 
u(x, ¥) = \2 eT ga | a ae 
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Zusammenfassung 


Fiir einen aus 7 Schichten verschiedener Warmeleitfahigkeit bestehenden un- 
endlichen Halbraum wird die Temperaturverteilung berechnet, die sich einstellt, 
wenn an der freien Oberflache ein beliebiger Temperaturverlauf vorgeschrieben 
wird. Eine geschlossene Darstellung ist im allgemeinen nicht méglich, doch gibt es 
Falle, wo die Lésung verhaltnismassig einfach wird. 


(Received: Juli 7, 1960.) 
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Numerische Methoden in der Meteorologie 


Von WALTER KuuHwn}), Ziirich 


Zusammenfassender Bericht 


1. Problemstellung, Vorlaufer 


Die Atmosphare empfangt Strahlungsenergie von der Sonne und von der 
Erdoberflache ; sie sendet selber Strahlung aus. An der fliissigen und festen Erd- 
oberflache werden zwischen Erde und Atmosphére Warme und Wasser ausge- 
tauscht. Die unterste Luftschicht unterliegt je nach ihrer Bewegung und je nach 
dem Zustand der Erdoberflache veranderlichen Reibungskraften. Abgesehen von 
diesen dusseren Einwirkungen, die sich quantitativ erfassen lassen, kann die 
Atmosphare als ein geschlossenes mechanisch-thermodynamisches System be- 
trachtet werden. Auf jedes Luftteilchen wirkt die Gravitation und die von der 
Rotation der Erde herriihrende Corioliskraft. 

Die Gesetze, welche die Zustandsanderungen im Innern der Atmosphare be- 
herrschen, sind im wesentlichen bekannt: 1. Die hydrodynamische Bewegungs- 


1) Schweizerische Meteorologische Zentralanstalt. 
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gleichung, aufspaltbar in drei Komponenten, 2. die Kontinuitatsgleichung, 3. die 
Zustandsgleichung idealer Gase, 4. der erste und zweite Hauptsatz der Thermo- 
dynamik, 5. die Gesetze der Phasenumwandlungen des Wassers, 6. die Gesetze 
des turbulenten Austausches von Impuls, Warme und Wasserdampf, 7. die Strah- 
lungsgesetze. 

Im Sinne der deterministischen Physik bestimmt ein momentaner Zustand 
alle spateren Zustande. Es stellt sich deshalb die Frage: Wie lassen sich bei hin- 
reichend genauer Kenntnis eines Anfangszustandes spatere Zustande berechnen ? 

BJERKNEs [1]?), der Schépfer der dynamischen Meteorologie, hat dieses Pro- 
blem bereits 1904 exakt formuliert und seine Schwierigkeiten klar erkannt. Sie 
liegen auf vier Gebieten: 

a) Liickenhafte und ungenaue Beobachtungen iiber den dreidimensionalen 

Anfangszustand, b) theoretisch-physikalische Schwierigkeiten: Erfassung der 
meteorologisch bedeutsamen Vorgange und Herausfilterung kleinrdumiger St6- 
rungen durch gezielte Vereinfachung der Grundgleichungen, c) mathematische: 
Integration eines komplizierten Systems partieller Differentialgleichungen mit 
unsicheren Randbedingungen, d) Arbeitsaufwand im Hinblick auf routinemassige 
Anwendung. 

: BJERKNES war sich klar dartiber, dass eine strenge analytische Lésung des 
Problems in voller Allgemeinheit kaum je in Frage kommen diirfte. Er empfahl, 
sich zuerst auf die rein hydrodynamische Seite zu beschranken und die Verquik- 
‘kung mit thermodynamischen Prozessen einer spaéteren Erweiterung vorzube- 
halten. Auch erkannte er, dass nur numerische oder graphische Integrationsmetho- 
den zum Ziele fiihren wiirden. BJERKNES wurde in der Folge durch seine Polar- 
frontlehre so stark in Anspruch genommen, dass er sich selber dem Problem der 
rechnerischen Wettervorhersage nicht mehr widmen konnte. 

Einen originellen und scharfsinnig angelegten Versuch unternahm in den 
Jahren 1911-1921 der Englander F.L. RicHarpson [2]. Leider stimmte das Er- 
-gebnis seiner miihevollen Rechnungen schlecht mit der Wirklichkeit tberein, 
ohne dass er und seine Zeitgenossen die Ursachen des Misserfolges v6llig klar er- 
kennen konnten. 

Es war naheliegend, dass man neben den hydrodynamischen Bewegungs- 
gleichungen vorab die Kontinuitétsgleichung heranzog, weil sie Anderungen 
der Dichte mit Vergenzen (Divergenzen und Konvergenzen) des Stromungsfeldes 
in Verbindung bringt. Da man iiber die vertikale Komponente der Luftbewegung 
a priori nichts wusste, versuchte man den Vergenzeffekt langs einer vertikalen 
Luftsaule vom Boden bis zur oberen Grenze der Atmosphare zu integrieren (so- 
genannte Tendenzgleichung), um die Druckaénderung am Boden zu berechnen. 
‘Dabei machte man aber zwei entmutigende Feststellungen: 1. Die Genauigkeit 
der Windmessungen reichte nicht entfernt zur Berechnung der Vergenzen aus. 
2. In einer vertikalen Luftsdule iiberlagern sich jederzeit gegensinnige Vergenzen, 
so dass der integrierte Effekt viel kleiner ist als die Beitrage einzelner Schichten. 

Damit war das Problem der Vorausberechnung von Wetterkarten einstweilen 
in eine hoffnungslose Sackgasse geraten, aus der es erst 1948 durch die bahn- 
brechenden Untersuchungen von CHARNEY und seinen Mitarbeitern befreit wurde. 

Inzwischen hatte die Aerologie, dh. die direkte Erforschung der freien Atmos- 
phare, durch die Schaffung eines weltweiten Netzes taglicher Radiosondierungen 
einen ungeahnten Aufschwung genommen. Aber auch die Dynamik der Luft- 
stromungen hatte einige Fortschritte aufzuweisen: 


2) Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis; Ziffern in runden 
Klammern auf die Formeln am Ende dieses Berichtes. 
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RossBy [3] hatte 1939 eine bestimmte.Wirbelgrosse, die sogenannte absolute 
Vorticity, in die Meteorologie eingefiihrt und mit Erfolg auf die Verlagerung langer 
Wellen in einer zonalen Grundstromung angewendet. Die absolute Vorticity ist 
definiert als Vertikalkomponente des absoluten Rotors, wobei wir unter dem ab- 
soluten Rotor den Drehvektor in bezug auf ein nicht mit der Erde rotierendes 
Inertialsystem verstehen. Demnach setzt sich die absolute Vorticity zusammen 
aus der relativen Vorticity, welche von der Drehung und Scheerung des horizon- 
talen Windfeldes herriihrt, und der Projektion des Erdrotations-Vektors auf die 
Vertikale; diese Komponente, Coriolis-Parameter genannt, hangt nur von der 
geographischen Breite ab. Fiir die absolute Vorticity gelten gewisse Erhaltungs- 
sitze, welche als Verallgemeinerungen des Helmholtzschen Wirbelsatzes fir 
kompressible Medien anzusehen sind. 

SUTCLIFFE [6] hatte 1947 mit Hilfe der Vorticity und der temperaturbedingten 
vertikalen Abstainde isobarer Flachen eine Theorie der Entwicklung thermisch 
asymmetrischer Zyklonen geschaffen. 


2. Das barotrope Modell als Grundlage 


Als die elektronischen Rechenautomaten aufkamen, legte der Mathematiker 
JouHN von NEUMANN 1946 einer Gruppe von Forschern am Institute for Advanced 
Study in Princeton die Frage vor, ob diese Gerate nicht neue Perspektiven fir die 
Vorausberechnung atmospharischer Stromungen erdffneten. Unter der Fiihrung 
von CHARNEY wandte sich diese Arbeitsgruppe dem Studium grossraumiger 
Str6mungen zu. CHARNEY [7] unterzog insbesondere die Grdéssenordnung aller in 
den dynamischen Gleichungen auftretenden Glieder einer eingehenden Unter- 
suchung, wobei ihre Abhangigkeit vom Ma8stab der betrachteten Phanomene 
klar zutage trat. 

Es ergab sich, dass die bei Schall- und Gravitationswellen wesentlichen Ver- 
genzen fiir grossradumige meteorologische Stromungen nicht charakteristisch sind. 
Diese schienen praktisch vergenzfret abzulaufen, und es musste also méglich sein, 
durch Unterdriickung der Vergenzen in den dynamischen Gleichungen die kleinen 
Storungen, den «meteorologischen Larm», herauszufiltern. 

Anderseits erwies sich die absolute Vorticity als eine in erster Naherung in- 
variante, fiir grossraumige Str6mungen in hohem Grade charakteristische Grosse. 
Auch spielten Temperaturunterschiede langs der Stromlinien dynamisch eine viel 
kleinere Rolle, als friiher angenommen worden war. In dem fiir meteorologische 
Zwecke interessanten Bereich des Spektrums atmospharischer Bewegungen durfte 
man also auch in erster Naherung Barotropie voraussetzen; darunter versteht man 
das Zusammenfallen der Flachen konstanter Temperatur mit Flachen konstanten 
Druckes. 

Wenn man kleinraumige Stérungen durch entsprechende Glattung des Strom- 
feldes eliminierte, schien sich die Luft in ihrem dynamischen Verhalten einer 
inkompressiblen Fliissigkeit zu nahern. Dieses vergenzfreie und quasi-barotrope 
Verhalten erwies sich vorziiglich als eine Eigenschaft der mittleren Tvoposphare, 
wahrend in den unteren und oberen Stockwerken der Troposphare gréssere Ab- 
weichungen gefunden wurden. Deshalb konzentrierte man die Berechnungsver- 
suche naturgemass zuerst auf die im Wetterdienst ohnehin gelaufige 500-Millibar- 
Flache, die im Durchschnitt etwa in 5,5 km Hoéhe iiber dem Meeresniveau liegt. 

Hier ist eine Bemerkung iiber das in der theoretischen Meteorologie iibliche 
Koordinatensystem am Platze. Als vertikale Koordinate verwendet man an Stelle 
der Hohe oder des Geopotentials mit Vorteil den Luftdruck bp, wahrend die Ebe- 
nen * = const und y = const wie in der Mechanik orthogonale Vertikalebenen 


- 
¢ 
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_ sind. Die Hohe z oder besser das Geopotential ® erscheint dann als abhangige 
_ Variable und wird meist dazu beniitzt, die Topographie einer isobaren Flache 
_ = const zu charakterisieren. Die Neigung isobarer Flachen ist iibrigens von der 
_ Grossenordnung 1:10000 und betragt nur in extremen Fallen 1:1000. Das (, y, 
_ p)-System ermoglicht die Elimination der Luftdichte aus samtlichen dynamischen 
_ Gleichungen [8]. 
Bei Untersuchungen, die sich auf die Erdatmosphare als Ganzes beziehen, 
‘muss man die horizontalen Koordinaten (x, y) durch sphdrische ersetzen; erstrek- 
ken sich die zugrundeliegenden Wetterkarten auf eine Hemisphare oder auf 
_ grosse Teile einer solchen, verwendet man in der Regel eine stereographische Pro- 
jektion und transformiert die dynamischen Gleichungen entsprechend. 
CHARNEY und Mitarbeiter zeigten nun auch, wie man die Vorticity prognos- 
tisch verwerten kann [9, 10]. Durch kreuzweises Differenzieren der beiden Be- 
wegungsgleichungen (1) nach y bzw. ¥ und Subtraktion, also durch Rotorbildung, 
erhalt man unter Beriicksichtigung der Kontinuitatsgleichung (2) die allgemeine 
Vorticitygleichung (3). Es sei darauf hingewiesen, dass hierin nur die Reibungs- 
bzw. Austauschkrafte und die von der vertikalen Windkomponente herriihrende 
_ Komponente der Corioliskraft vernachlassigt wurden; die letztgenannte Vernach- 
ldssigung ist bei grossraumigen Stro6mungen belanglos, da die Vertikalgeschwin- 
digkeit im grossrdumigen Durchschnitt rund 1000 mal kleiner ist als die horizon- 
tale. 
Im vergenzfreien, barotropen Fall verschwindet die rechte Seite von (3) und 
das Vorticitygesetz reduziert sich auf einen einfachen Erhaltungssatz (4). Er be- 
sagt, dass die absolute Vorticity einer Luftmasse im Laufe der Bewegung erhalten 
bleibt, d.h. sich mit dem Wind verlagert. Dabei ist zu beachten, dass infolge der 
_Breitenabhangigkeit des Coriolisparameters die relative Vorticity im allgemeinen 
nicht konstant bleibt: wandert eine Luftmasse gegen den Aquator, so nimmt ihre 
Drehbewegung relativ zur Erde in zyklonalem Sinne zu, und umgekehrt. 
Nun ergibt sich aber eine weitere Vereinfachung aus der unter gewissen Voraus- 
setzungen zuldssigen geostvophischen Naherung. Als geostrophischen Wind be- 
-zeichnet man einen Wind, der sich bei stationadrer Str6mung, d.h. bei Abwesenheit 
von Reibung und Beschleunigungen, aus dem Gleichgewicht von horizontaler 
Druckkraft und Corioliskraft ergibt. Aus den aerologischen Karten weiss man seit 
langem, dass der Wind in der freien Atmosphare, und zwar namentlich in der 
‘Mitte der Troposphare, im grossen und ganzen dem geostrophischen Naherungs- 
gesetz sehr gut gehorcht. Merkliche Abweichungen treten nur bei scharfer Kriim- 
mung, Diffluenz oder Konfluenz der Stromlinien auf. 

Der geostrophische Wind weht parallel zu den Isohypsen (Héhenkurven) der 
isobaren Flachen mit einer Starke, die nur von der Neigung der isobaren Flachen 
und von der geographischen Breite abhangt. 

Bei Verwendung der geostrophischen Naherung kénnen wir die oft nur un- 
genau gemessenen Winde durch einen horizontalen Vektor ersetzen, der sich in 
einfacher Weise aus der viel genauer bestimmbaren Topographie isobarer Flachen 
berechnen lasst. So ergibt sich aus (4) eine recht einfache Prognosengleichung fiir 
die zeitliche Anderung des Geopotentials (6); sie hat die Form einer Poissonschen 
Differentialgleichung, wobei die rechte Seite ausschliesslich von topographischen 
Eigenschaften (Kriimmung und Neigung) der betrachteten isobaren Flache und 
von der geographischen Breite abhangt. Die Auflésung dieser Gleichung nach 
0@/dt hat eine gewisse Verwandtschaft mit der Berechnung des elektrostatischen 
Feldes aus einer gegebenen Raumladung, ist aber hier nur fiir ein zweidimensiona- 
les Problem durchzufiihren. 
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Ist 0/ot einmal fiir jeden Punkt der Flache berechnet, so geniigt Multiplika- 
tion mit einem bestimmten Zeitintervall At und Addition zum Ausgangswert fir 
die Berechnung des Geopotentials am Ende des Zeitintervalles (7). Der Zeit- 
schritt At ist allerdings begrenzt durch die Forderung, dass die rechte Seite der 
Prognosengleichung (6) wahrend dieser Zeitspanne innerhalb der Fehlergrenzen 
konstant bleibe. Bei den in Betracht kommenden Feldanderungen betragt die 
héchstzulassige Prognosendauer fiir einen Schritt 2 Stunden. Eine Prognose aut 
24 Stunden bedingt deshalb 12 malige Wiederholung des Verfahrens, wobei die 
bei jedem Schritt erhaltene ®-Verteilung als Ausgangsverteilung fiir den nachsten 
Schritt zu verwenden ist. 

Fiir die Durchfiihrung solcher Aufgaben stehen sowohl numerische wie gra- 
phische Verfahren zur Verfiigung. 

Die hauptsachlich von FjGrtTorr [12] und EsroguE [16] entwickelten gra- 
phischen Verfahren sind den numerischen im Prinzip ahnlich; durch einen beson- 
deren Kunstgriff konnte der Zeitschritt fiir eine graphische Prognose ohne starke 
Beeintrachtigung der Genauigkeit sogar auf 24 Stunden erhéht werden; auch ist 
der finanzielle Aufwand fiir graphische Prognosen bedeutend geringer als der- 
jenige fiir numerische Prognosen mit Hilfe elektronischer Rechenanlagen. Da die 
graphischen Verfahren aber das Zeichnen einer grossen Zahl von Hilfskarten er- 
fordern, haben sie sich fiir Routinezwecke neben den numerischen nicht durch- 
setzen kénnen. Im folgenden befassen wir uns ausschliesslich mit den numerischen 
Methoden. 


3. Numerische Berechnung 


Um irgend ein physikalisches Feld numerischer Behandlung zuganglich zu 
machen, muss man seine Zahlenwerte in regelmassig angeordneten diskreten 
Punkten des Raumes kennen; das heisst, man muss den kontinuierlichen Raum 
durch ein sogenanntes Gittey unterteilen. Die Wahl der Gitterkonstanten oder 
Maschenweite richtet sich nach der rdumlichen Ausdehnung der zu untersuchen- 
den Feldschwankungen. Fiir die Bearbeitung grossrdumiger atmospharischer 
«Stérungen» (Wirbel und Wellen) hat sich eine horizontale Maschenweite von etwa 
400 km als zweckmassig erwiesen [25]. 

Nun sind aber die aerologischen Beobachtungsdaten zunachst nicht in den 
Punkten eines regularen Gitters, sondern an unregelmdssig verteilten Stationen 
gegeben. (Eine Anzahl stationarer Wetterschiffe sorgt fiir Sondenaufstiege von den 
Ozeanen aus; tiber den Landgebieten ist die Stationsdichte wesentlich grésser). 
Jeder Prognose hat deshalb eine Analyse des Beobachtungsmateriales voranzu- 
gehen, durch welche man die wahrscheinlichen Werte der meteorologischen Va- 
riablen in den Gitterpunkten fiir einen bestimmten Ausgangstermin festlegt. In 
der konventionellen Meteorologie erfolgt diese Analyse, indem man in die Karte 
durch Interpolation zwischen den vorliegenden Stationswerten von Hand még- 
lichst glatte Isolinien zeichnet. Beim Zeichnen der Isohypsen einer isobaren Flache 
nimmt man auf das geostrophische Naherungsgesetz Riicksicht, indem man ver- 
sucht, die Tangentenrichtung und den Abstand der Kurven méglichst gut in Uber- 
einstimmung mit den gemessenen Winden zu bringen. Offensichtlich falsche 
Meldungen sind dabei nach Méglichkeit zu korrigieren oder zu ignorieren. Nach 
Beendigung einer solchen Analyse sind fiir numerische Zwecke die Werte in den 
Gitterpunkten abzulesen, das heisst durch Schatzung zu interpolieren. 

Bei diesem Arbeitsprozess ist eine gewisse persénliche Willkiir namentlich in 
Gebieten, wo die Beobachtungsstationen weit auseinanderliegen, nicht zu um- 
gehen. Um solche personliche Fehler auszuschalten, haben CRESSMAN und andere 
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‘ mote Analysenmethoden [24] erfunden, welche mit Hilfe elektronischer Re- 
_chenautomaten nach dem Prinzip der kleinsten Quadrate und mit geostrophi- 
_ scher Beriicksichtigung des Windes auf rein maschinellem Wege zu den Ausgangs- 
_ werten fiihren. Diese Ausgangswerte werden dann der Prognosenmaschine «ver- 
_ fittert». Im einfachsten Fall der barotropen Prognose einer 500-Millibar-Karte 
_ besteht das Ausgangsmaterial aus der zweidimensionalen @-Verteilung, dh. aus 
_ den Werten des Geopotentials der 500-Millibarflache an einem bestimmten Termin. 
' Die im Speicherwerk des Computers ein fiir allemal festgelegte Rechenvor- 
_ schrift stiitzt sich auf die im Anhang mitgeteilten Formeln. Dabei sind Differen- 
_ tialquotienten durch entsprechende Differenzenquotienten anzunahern. Die Ma- 
_ schine hat also im wesentlichen die ersten und zweiten Differenzen des Geopoten- 
; tials ® nach den Koordinatenrichtungen ¥# und y und daraus die relative Vorticity 
_ ¢ zu bilden; durch Addition des zeitlich konstanten Coriolisparameters / wird die 
_ absolute Vorticity 4 erhalten. Von ihr sind wiederum die ersten Differenzen zu 
_ bilden, damit der Jacobi-Ausdruck J(7, ®) berechnet werden kann. Fiir die Auf- 
_ lésung der Poisson-Gleichung (6) nach 0@/d¢ dient ein Relaxationsverfahren [4, 11, 
_ 13]; der Rest besteht aus elementaren Operationen gemass (7). 
Bisher haben wir die Randwertprobleme iibergangen. Bei jeder Differenzbil- 
dung wiirde das nutzbare Gitter an allen Randern um eine Maschenbreite ver- 
_ kleinert, wodurch das Prognosengebiet sehr rasch zusammenschrumpfen wiirde. 
_ Man ist also gezwungen, an den Randern bei jeder Zwischenoperation willkiirliche 
_ @-Werte vorzugeben. Dies ist namentlich am Westrand eines Kartengebietes eine 
sehr bedenkliche Massnahme, da durch die vorherrschenden Westwinde der Ein- 
_fluss der Randwerte rasch ins Innere des Kartengebietes fortgepflanzt wird. Um 
diese Fehler zu vermeiden, gibt es nur eine radikale Loésung: Das Kartengebiet 
muss womdglich auf eine von Parallelkreisen begrenzte, um die ganze Erde herum- 
laufende Zone oder noch besser auf eine Kugelkalotte, ja Hemisphare ausgedehnt 
_ werden. Da in den Tropen @ erfahrungsgemass nur wenig andert, ist dort die Vor- 
gabe konstanter ®-Werte zulassig. 

Heute wird in den USA bei numerischen Routineprognosen ein auf eine 
stereographische Projektion gelegtes, achteckig begrenztes Gitter verwendet, das 
fast die ganze Nordhalbkugel bis zum 10. nordlichen Breitengrad bedeckt. Dieses 
im Innern in lauter Quadrate unterteilte Gitter weist nahezu 2000 innere Netz- 
punkte auf. 

Das Ergebnis einer numerischen Prognose besteht aus Zahlen, die von der 
Maschine in richtiger Anordnung auf eine stereographische Karte gedruckt wer- 
den, worauf ein Zeichner bloss noch Isolinien des neuen (vorausberechneten) 
Geopotentials (=: Stromlinien) einzuzeichnen hat. Diese Prognosenkarte wird 
dann durch drahtlose Facsimile-(Bildfunk-) Ubertragung in kiirzester Zeit allen 
Interessenten zugestellt. Der ganze Arbeitsprozess von der Eingabe der Stations- 
beobachtungen bis zur Ausgabe der Prognosenkarte dauert etwa eine Stunde und 
lauft praktisch ohne menschliches Dazutun ab. 


4. Barokline und ageostrophische Modelle 


Es ist klar, dass die barotropen Gleichungen nicht streng richtige Vorhersagen 
liefern kénnen; denn sie verlangen ja Konstanz der Wirbelintensitat, wahrend 
anhand von Wetterkarten immer wieder festgestellt werden kann, wie sich Wirbel 
von einem Tag zum nachsten verstairken, abschwachen, auflésen oder neu bilden. 
Nach der Entwicklungstheorie von SuTcLIFFE [6] liegt die Ursache solcher Wir- 
belanderungen wenigstens zum Teil in baroklinen Effekten, d.h. in isobaren Tem- 


peraturunterschieden. 
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Um ihnen Rechnung zu tragen, muss man das erste Glied auf der rechten Seite 
der allgemeinen Vorticitygleichung (3) beriicksichtigen. Ein Zusammenhang — 
zwischen Temperaturanderungen und Vertikalgeschwindigkeit wird durch den | 
ersten Hauptsatz der Warmelehre hergestellt, wobei man sich anfanglich auf 
adiabatische Vorgange beschranken darf. 

Zur Erfassung isobarer Temperaturunterschiede geniigt eine einzige isobare— 
Flache als Ausgangsmaterial nicht mehr; man braucht mindestens zwei solcher 
Flachen, wobei ihr vertikaler Abstand iiber die Mitteltemperatur der dazwischen 
liegenden Luftschicht Aufschluss gibt. Bei Verwendung von zwei isobaren Fla-_ 
chen kann man die Temperatur in der Horizontalen varieren lassen, muss aber | 
einen einheitlichen vertikalen Temperaturgradienten voraussetzen. Will man auch 
diesen von Ort zu Ort varieren lassen, muss eine dritte isobare Flache herangezo- 
gen werden; eine vierte erlaubt Anderungen des vertikalen Temperaturgradien- 
ten auch in vertikaler Richtung, usw. So werden der theoretischen Modellatmos- 
phare sukzessive weitere Freiheitsgrade eingerdumt, wodurch aber das zu verar- 
beitende Material immer umfangreicher wird [13, 14]. 

Bei diesen Mehrschichtmodellen liefert die Maschine als Ergebnis die vorausbe- 
rechnete Topographie samtlicher isobaren Flachen, deren Anfangstopographie ihr > 
mitgeteilt wurde. Die Anforderungen an die Speicherkapazitat der Computers 
sind dadurch gegeniiber der barotropen Rechnungsweise vervielfacht worden. 

Interessanterweise liefern die baroklinen Modelle fiir die Topographie der 500- 
Millibar-Flache im allgemeinen kaum bessere Prognosen als das barotrope Mo- 
dell — ein indirekter Beweis fiir das quasi-barotrope Verhalten der mittleren Tro- 
pospharenschichten. Fiir héher oder tiefer gelegene Schichten ist dagegen die 
barokline Berechnung unumganglich. 

Wie primitiv das barotrope Modell ist, geht schon aus der Tatsache hervor, 
dass bei barotroper Temperaturschichtung der Wind von der Héhe unabhangig 
ware. In Wirklichkeit sieht aber das Stromfeld in verschiedenen Héhen oft ganz 
anders aus. 

Alle bisher besprochenen Modelle leiden aber an einer gemeinsamen Schwache, 
namlich an der geostvophischen Naherung. Diese ist zwar sehr bequem, aber physi- 
kalisch durchaus nicht gerechtfertigt. Ein streng geostrophischer, also parallel zu 
den Niveaulinien der isobaren Flachen wehender Wind bewirkt naémlich keine 
Dichteverlagerung und kénnte infolgedessen iiberhaupt zu keinen Anderungen 
der Druckverteilung Anlass geben. Dass man mit den geostrophischen Modellen 
Anderungen des Druck- und Windfeldes vorhersagen kann, liegt nur an gewissen 
Inkonsequenzen in der Anwendung der geostrophischen Naherung. 

Es bedeutet einen prinzipiellen Fortschritt, dass es 1955 CHARNEY [15], 
Bottn [20] und anderen [19] gelang, die geostrophische Fessel durch die soge- 
nannte Balance-Equation (8) zu sprengen. 

Bei der Ableitung dieser Gleichung geht man von der Tatsache aus, dass gross- 
rdumige StrOmungen in guter Naherung vergenzfrei sind, sich also durch eine 
Stromfunktion Y darstellen lassen. Durch Einsetzen der Stromfunktion in die 
urspriinglichen hydrodynamischen Bewegungsgleichungen (1) und Divergenz- 
bildung erhalt man nach leichter Umgruppierung die Balance-Equation (8). Sie 
stellt eine Differentialgleichung zur Bestimmung der zunachst unbekannten 
Stromfunktion YW dar, falls die Topographie ® einer isobaren Flache gegeben ist, 
und kann mittels eines Relavationsverfahrens [4] numerisch gelést werden. Die 
aus der Stromfunktion sich ergebenden Winde besitzen im allgemeinen eine 
kleine ageostrophische Komponente; die Stromlinien Y= const verlaufen nicht 
mehr genau entlang den Isohypsen. 
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__ Zur Beriicksichtigung ageostrophischer Effekte ist das Rechnungsprogramm 
der numerischen Wetterprognose in der Weise abzudndern, dass zunachst von ® 
_ausgehend die Stromfunktion Y berechnet wird; alle weiteren Operationen wer- 
den dann mit dieser Stromfunktion ausgefiihrt; erst beim Schlussresultat geht 
man wieder zur Topographie iiber. Dieses Verfahren hat sich namentlich in den 
_ Fallen bewahrt, wo mit grossen Beschleunigungen zu rechnen ist, zum Beispiel bei 
_tropischen Wirbelstiirmen [21]. 

: Weitere Verbesserungen sind zu erwarten von der Beriicksichtigung der Ge- 
_lande-, Reibungs- und Strahlungseinfliisse sowie der Kondensationsprozesse. 


a 


5. Bisherige Erfolge und Zukunftsaussichten 


s 


; Unsere Darstellung nimmt im wesentlichen Bezug auf die Arbeiten der Joint 
Numerical Weather Prediction Unit (JNWP) in Suitland (Maryland, USA). Dieses 
von der US-Air-Force, der amerikanischen Marine und dem staatlichen Wetter- 
bureau gemeinsam finanzierte und betriebene Institut macht seit 1955 tagliche 
_Routineprognosen; es verwendet einen «Digital Computer» vom Typ IBM-704, 
der demnachst durch eine Anlage vom Typ IBM-7090 ersetzt wird, also durch 
eine der gréssten bisher gebauten Rechenmaschinen. 
Die Produkte der JNWP werden laufend durch Facsimile-Funkiibertragung 
-ausgestrahlt. Das gegenwartige Tagesprogramm umfasst: Eine 72-stiindige Voraus- 
berechnung der 500-Millibar-Karte (Niveau 51% km), eine 24-stiindige und zwei 
-in Abstanden von 12 Stunden ausgegebene 36-stiindige Prognosen der 700-Milli- 
_bar-Karte (Niveau 3 km), je zwei 12-, 24- und 36-stiindige Prognosen der 300- 
-Millibar-Karte (Niveau 9 km), eine 24-stiindige Prognosenkarte der grossraumigen 
'Vertikalbewegung. 
Die 300-Millibar-Prognosen sind von unmittelbarer Wichtigkeit fiir den Diisen- 
Luftverkehr. Vertikalbewegungskarten sollen eine Abschatzung der Bewolkung 
und Niederschlage erméglichen. 

Auf Grund von Prognosen der Vertikalbewegung und der im Ausgangszeit- 
punkt bekannten dreidimensionalen Feuchtigkeitsverteilung ist es tibrigens be- 
reits gelungen, quantitative Niederschlagsprognosen zu machen, wobei sogar die 
durch das Gelande bedingten Aufwinde naherungsweise in Rechnung gestellt wer- 
den konnten. Diese Niederschlagsprognosen gaben die grossraumige Verteilung 
des Niederschlages im allgemeinen recht zutreffend wieder, wogegen Instabilitats- 
erscheinungen regionalen und lokalen Charakters (Schauer und Gewitterregen) 
natiirlich nicht erfasst werden konnten. 

Es gibt auch schon brauchbare Verfahren zur Voraussage des Meereszustandes. 
‘Hierbei wird der Wellengang auf Grund von numerischen Vorhersagen des Windes 
in der untersten Luftschicht berechnet. 

Diese Beispiele zeigen, dass neuerdings numerische Methoden auch auf andere 
Gréssen als den Wind oder das Geopotential isobarer Flachen angewendet werden. 
Im Ubrigen sei jedoch hervorgehoben, dass numerische Prognosen vorlaufig nur 
ein grossydumiges Bild iiber die Verteilung des Windes oder allenfalls eines anderen 
meteorologischen Elementes geben. Die Interpretation dieses Bildes im Hinblick 
auf vegionale Wettervorhersagen bleibt nach wie vor Aufgabe erfahrener Meteoro- 
logen, die ihre Regeln teils aus physikalischen Gesetzen, teils aus dem Studium 
vergangener Wetterabldufe herleiten. Wer aber im Wetterdienst arbeitet, weiss, 
wieviel fiir die regionale Wettervorhersage schon mit einer zuverlassigen Prognose 
des grossrdumigen Strémungsfeldes gewonnen ist. 

Wie steht es heute mit der Giite numerischer Prognosen ? Um diese objektiv zu 
priifen, hat man Korrelationskoeffizienten zwischen vorausberechneten und ein- 
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getroffenen Anderungen bestimmter Elemente an den Punkten eines Gitters 


berechnet. Dabei ergaben sich bisher folgende Resultate: : 
Im 500-Millibar-Niveau sind die numerischen Eintages-Prognosen den empi- 
rischen ebenbiirtig; bei langerer Prognosenfrist (2 oder 3 Tage) sind die numeri- 


schen Prognosen den empirischen deutlich iiberlegen. In anderen Stockwerken ; 


der Atmosphare, insbesondere am Erdboden, ist die Zuverlassigkeit der empiri- 


schen Methoden von den numerischen bisher erst knapp erreicht worden. Doch | 


besteht kaum ein Zweifel, dass mit der weiteren Verfeinerung der Modelle die 
numerischen Methoden mehr und mehr den Vorrang erhalten werden. 

Ihre kiinftige Entwicklung richtet sich aber noch auf ein anderes Anwendungs- 
gebiet: Die Mittel- und Langfristprognosen. 

Unter «Mittelfrist» verstehen wir in der Meteorologie einen Zeitraum, der mit 
der Dissipationszeit einer normalen Depression vergleichbar ist (etwa 5—7 Tage). 
Fiir solche Zeitspannen ist eine exakte dynamische Vorhersage auf Grund des 
Anfangszustandes heute nicht méglich. Dagegen scheint es auf Grund von Be- 
ziehungen zwischen raéumlichen und zeitlichen Mitteln, die u.a. von 6sterreichi- 
schen Meteorologen [18] untersucht wurden, zu gelingen, generelle Aussagen tiber 
die wesentlichen Ziige des atmospharischen Strémungsbildes am Ende eines 
mittelfristigen Zeitraumes zu machen [22]. 

Unter Langfristprognosen verstehen wir Voraussagen auf Monate, Jahreszei- 
ten oder Jahre hinaus. Bei solchen Zeitraumen spielen Einzelheiten des Anfangs- 
zustandes, soweit man sie heute erfassen kann, keine Rolle mehr. Es handelt sich 
vielmehr darum, die Abweichungen der sogenannten <«allgemeinen Zirkulation» der 
Erdatmosphare vom mittleren Verlauf vorauszuberechnen. 

Nach heutiger dynamischer Auffassung ist die allgemeine Zirkulation das, was 
von einem momentanen Stro6mungszustand der Erdatmosphare iibrig bleibt, wenn 
man durch eine starke Glattung alle Stdrungen von der Gréssenordnung kleiner 
bis mittlerer Depressionen weggeschliffen hat. Sie zeigt zyklonale und antizyklo- 
nale Wirbel oder auch Wellen grossen Ausmasses, welche sich einer zonalen Grund- 
stro6mung tiberlagern und fiir den meridionalen Warmeaustausch zwischen den 
verschiedenen Breitenzonen verantwortlich sind. Bei Mittelbildung iiber langere 
Zeitraume (Jahre) bleibt abgesehen von Land-Meer-Effekten nur die zonale Grund- 
stro6mung erhalten, die das dynamische Geschehen nur zu einem kleinen Teil wider- 
spiegelt. Die Theorie der allgemeinen Zirkulation steckt noch in ihren Anfangen. 

Die numerischen Methoden sind berufen, Licht in diese Theorie zu bringen ; 
denn sie geben uns ein Mittel in die Hand, mit theoretischen Modellen zu experi- 
mentieven [17], was in natura nicht méglich ist. Wir kénnen der Rechnungsma- 
schine z.B, gewisse Grundtypen einer allgemeinen Zirkulation eingeben und 
untersuchen, wie sich die Zirkulation im Laufe der Zeit unter bekannten Ausseren 
Einwirkungen verdndert. Dabei miissen selbstverstandlich Strahlung, Reibung, 
Kondensationsprozesse und topographische Einfliisse in Rechnung gestellt werden; 
vielleicht sind aber auch extraterrestrische Einfliisse wie die Schwankungen der 
Sonnenstrahlung zu beriicksichtigen. 

Verheissungsvolle Ansatze zur Lésung dieser Probleme sind bereits vorhanden. 
Die Ergebnisse sind zum Teil iiberraschend. So scheint die iiber halbe Erdqua- 
dranten sich erstreckende «Grossturbulenz» mit ihren thermisch asymmetrischen 
Wirbeln energetisch ganz anders zu wirken als die bekannte Kleinturbulenz. 

Parallel zu diesen numerischen Versuchen werden in Laboratorien phystha- 
lische Modellversuche mit rotierenden, beheizten Fliissigkeiten gemacht, wobei sich 
zum Teil verbliiffende Ubereinstimmungen mit den aus Wetterkarten bekannten 
Stro6mungsbildern und auch mit den numerischen Experimenten ergaben [23]. 
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Fir die Durchfiihrung von Routineprognosen geniigen beim heutigen Stand 
der Ubermittlungstechnik einige wenige Zentren, zum Beispiel eines pro Konti- 
nent. Zur Entwicklung der theoretischen Grundlagen sind dagegen Forscher aller 


_ Lander aufgeru fen. 


Die in unserer Ubersicht besprochenen Methoden verdanken wir in der Haupt- 


_ sache nordamerikanischen Universitaten und F6rschungsinstituten; wesentliche 


Beitrage stammen von schwedischen, norwegischen, englischen, deutschen und 
japanischen Forschern. Russland geht mit den Methoden von I.A. KiBEL [5] 


_ eigene Wege. In manchen Landern wird auch versucht, die Resultate numerischer 


Berechnungen durch empirische Korrekturen zu verbessern. Das Entscheidende 


‘sind aber wohl die Entdeckungen von CHARNEY [7, 9] und seinen Mitarbeitern, 


durch welche eine physikalische Grundlage fiir objektive, mathematische Pro- 
gnosen des atmospharischen Stro6mungsfeldes auf Grund hydrodynamischer Ge- 


_setze geschaffen wurde. 


Symbole 
yy rechtwinklige, horizontale Koordinaten; 
p Luftdruck; 
uU, U Geschwindigkeitskomponenten nach x, y (Horizontalwind) ; 
Po = ie sog. generalisierte Vertikalgeschwindigkeit ; 
® = © (p; x, y; t) Geopotential; 
z = e +u ¢ +0 +o : individuelle (das heisst massengebundene) 
dt Ot Ox oy 0 


Ableitung nach der Zeit; 


V= (= ¥) Nabla-Operator auf einer isobaren Flache; 
VIP 


OX -a00 
02 02 5 * oo 4 
a ( 73 JC ay i Laplace-Operator auf einer isobaren Flache; 


da 0b da db 
INC ale dy Ody Ox 


Jacobi-Operator; 


: Fah aed sane etl oe OV a OME 

W Stromfunktion, definiert durch 4 = aad ,. Jr Fan 
1 0@ 1 0@ ; aes 
Ug = — Foy”? Vg = a geostrophischer Wind; 


f= 22-sing Coriolis-Parameter ; 


Q Rotationsgeschwindigkeit der Erde; 
2 geographische Breite; 
(ie eee = y*¥ relative Vorticity; 

Ox oy 


n=¢+f absolute Vorticity. 
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Wichtigste Formeln . 
Ou Ou Ou Ou 0® | 


+ u 0 +o stags, oS +fU1] Horizontalkomponenten der 
(1) a of ey a fs Eulerschen Bewegungsgleichung _ 
dv dv Ov ov,  0@ im (x, v, p)-System. 
= —fu 
‘Viele Sake daar) oe ee 


(2) ee ale - — _ = 0 Kontinuitatsgleichung im (4, y, p)-System. 


4 dn Ow du dw dv Ow 
Ue dice Soar. BPs Oy. OP uae 


Allgemeine Vorticitygleichung. 


ay ie os oD Barot Vorticitygleichung. 

(4) Fite 0 oder ai Cie, Uv ay o Op arotrope icityg g 

Ow On ; ; 
(5) ears = J(n, ¥)—o “OP. mit = 72% +f Barotrope Prognosengleichung. 
(6) ple S IGy QD) wile G= nk y?® + f Genaherte geostrophisch-barotrope 

of f Prognosengleichung. 

; 0® 
(7) We + At) = W(t) + At, oe bzw. @(¢ + At) = Ot) + At. 
GALE GEGE 0? \2 : 

(8) {VY+7fve+ 2 [ ar (5 a = V*® Balance-Equation. 
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Résumé 


Le but de cet article est d’exposer a des physiciens et des mathématiciens les 
problémes et principes de la prévision dynamique du temps et leur solution par 
des méthodes numériques. 


4 


Summary 


The purpose of the present article is to acquaint physicists and mathemati- 
cians with the problems and principles of dynamical weather prediction and their 
solution by numerical methods. 


(Eingegangen: 21. Juli 1960.) 
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Tagungen iiber Elektronenmikroskopie | 


1. Die Jahrestagung der Deutschen Gesellschaft fiir Elektronenmikroskopie é: Vw 
findet vom 24. bis 27. September 1961 in Kiel im Anatomischen Institut (Dir. Prof. — 
W. BarcMann) statt. Anfragen sind zu richten an: Dr. H. KEHLER, p. Adr. Farb- 
werke Hoechst AG., Frankfurt (Main)-Héchst. : 

2. Der alle vier Jahre stattfindende Internationale Kongress der International 
Federation of Electron Microscopy Societies wird auf Einladung der Amerikani- 
schen Gesellschaft fiir Elektronenmikroskopie vom 29. August bis 5. September 
1962 in Philadelphia stattfinden. Anfragen sind zu richten an: Fifth International 
Congress for Electronic Microscopy, 7701 Burholme Avenue, Philadelphia 11, 
Pennsylvania, USA. ; H. KEHLER 
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Games and Decision: Introduction and Critical Survey. Von R. D. LucE 
und H. Ratrra (John Wiley & Sons, Inc., New York 1957). 509 S., 59 Fig.; $8.75. 

Nach den Worten der Verfasser ist dies ein Buch iiber Spieltheorie und keine 
Darstellung der Theorie selbst. Dies kommt vor allem darin zum Ausdruck, dass 
es fast keine Beweise enthalt — auch keinen fiir die Existenz der Lésung des Zwei- 
Personen-Nullsummen-Spiels. Das Buch ist somit vor allem solchen Lesern zu 
empfehlen, fiir welche mathematische Beweise die Darstellung mehr verdunkeln 
als erhellen. — Nach einer allgemeinen Einleitung wird die Normalform von KUHN 
dargestellt, welche Spiele als Baume (im Sinne der Graphentheorie) interpretiert. 
Dann wird das Zwei-Personen-Spiel und das allgemeine v-Personen-Spiel behandelt, 
wofiir der Lésungsbegriff von v. NEUMANN-MORGENSTERN ausfihrlich erlautert 
und mit Experimenten verglichen wird, welche bei RAND gemacht worden sind. 
In weiteren Paragraphen wird das Problem einer individuellen Entscheidung bei 
unvollstandiger Information und jenes von Gruppenentscheidungen besprochen. 
Einige Anhange erweitern die Darstellung des Hauptteiles nach der mathemati- 
schen Seite. Es wird so das Minimaxtheorem exakt formuliert und sein Beweis 
auf den Beweis eines Fixpunktsatzes zuriickgefiihrt, und es werden die verschie- 
denen geometrischen Interpretationen des Zwei-Personen-Nullsummen-Spiels be- 
sprochen sowie auf den Zusammenhang mit der linearen Programmierung ein- 
gegangen. E, SPECKER 


Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von 
S. FLUGGE, Band 53: Astrophysik IV: Sternsysteme (Springer-Verlag, Berlin 1959). 
Mitarbeiter: F. K. Edmondson, B. Lindblad, Jj. HS Oort, HS Hoge, Ro AY Browse 
B. Y. Mills, G. de Vaucouleurs, F. Zwicky, J. Neyman, E. L. Scott, G. C. McVittie, 
O. Heckmann und E. Schiicking. 565 S., 189 Abb.; DM. 142.-. 

Der vorliegende letzte der die Ergebnisse der astrophysikalischen Forschung 
enthaltenden vier Bande behandelt die Sternsysteme und die Kosmologie. Das 
galaktische System nimmt 240 Seiten in Anspruch, die extragalaktischen Systeme 
200 und die Kosmologie 100 Seiten, Der Band gibt mehr eine lockere Aneinander- 
reihung ausgewahlter aktueller Kapitel als eine kompakte und geschlossene Dar- 
stellung. Diese Kapitel behandeln die Kinematik und Dynamik des galaktischen 
Systems, die radioastronomische Erforschung der Milchstrasse, die Sternhaufen, 
die «Punktquellen», die radioastronomische Erforschung, Klassifikation und Be- 


led 


Ps 
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_ schreibung der extragalaktischen Systeme, die Haufenbildung und die raéumliche 
_ Verteilung derselben, Entfernung- und Altersbestimmungen und schliesslich die 
_ kosmologischen Theorien. 
= Der Band enthalt eine reiche Fiille neuer und neuester Beobachtungsergebnisse 
und vorziigliches Bildmaterial. Die Darstellung geht je nach Gegenstand und Autor 
_ von der reinen Beschreibung, welche vorherrscht, bis zur reinen Theorie. Vielfach 
_ geben die Beitrage mehr eine Ubersicht der vom Verfasser und seinen Mitarbeitern 
_ gepflegten Forschungsrichtung als einen gesamthaften Uberblick. Einzelne Gebiete 
_ haben in neuester Zeit bereits eine monographische Darstellung, zum Teil von den- 
_ selben Autoren, erfahren, denen der Handbuchartikel kaum wesentlich Neues bei- 
 zufiigen hat. Andere aktuelle Gebiete, fiir welche dies nicht zutrifft, sind auch im 
vorliegenden Band iibergangen worden. Insbesondere vermisst man einen Beitrag 
uber interstellare Materie und magnetische Felder des galaktischen Systems. Die 
drei radioastronomischen Beitrage, welche zusammen nur 100 Seiten umfassen, 
hatte man gerne in grésserer Breite gesehen. So mag zwar der Leser nicht alles 
finden, was er sucht, aber was er findet, ist von bestem Standard und von der Héhe 
der Forschung geschrieben. Ein Vergleich mit dem vor einem Vierteljahrhundert 
erschienenen Handbuch der Astrophysik zeigt die erstaunliche Fiille neuer Kennt- 
_ nisse, denen leider keine entsprechende Mehrung des Verstehens folgte. Die Lésung 
des kosmologischen Problems ist noch in weiter Ferne, die meisten radioastronomi- 
_schen Beobachtungen sind unverstanden und die Entwicklung der Sternsysteme 
ratselhaft, wissen wir doch nicht einmal mit Sicherheit, in welcher Richtung sie sich 
drehen! M. WALDMEIER 


Figures of Equilibrium of Celestial Bodies. Von ZDENEK Kopat (The 
University of Wisconsin Press, Madison 1960). 153 S., 3 Fig.; $3.-. 

Die Theorie der Gleichgewichtsfigur einer kompressiblen Fliissigkeit unter der 
Wirkung der Eigengravitation wurde von CLAIRAUT, LEGENDRE und LAPLACE im 
wesentlichen bis zur 1. Ordnung des Einflusses der Deformation entwickelt. Im 
vorliegenden Buch gibt Kopar zundchst eine Zusammenfassung der Clairautschen 
Theorie, die explizite Ausdriicke fiir das Potential liefert und die 4ussere Form 
als Aquipotentialflache definiert. Im besonderen werden die Lésungen der Clai- 
rautschen Gleichung fiir spezielle Dichteverteilungen angegeben und ebenso die 
iterative Lésung im Falle beliebiger Dichteverteilungen, wie sie durch LjAPUNOV 
und Kopat entwickelt wurde, dargestellt. Das Hauptziel des Buches ist die erst- 
malige Darstellung der vollstandigen hydrostatischen Theorie 2. Ordnung. Es 
werden zu diesem Zwecke zunachst die expliziten Ausdriicke fiir die Aquipotential- 
flachen und das dussere Gravitationspotential abgeleitet fiir den Fall von Gezeiten- 
stérungen und Rotationseffekten bei beliebiger Dichteverteilung. Speziell werden 
dann das homogene Modell und das Rochesche Modell betrachtet. Ebenso werden 
die Tragheitsmomente berechnet. Im weitern werden die recht komplizierten 
Terme bestimmt, die aus der Wechselwirkung von Gezeiten und Rotation auftreten 
und hier nicht mehr vernachlassigt werden diirfen. Schliesslich untersucht KopaL 
noch die nichtradialen Schwingungen, insbesondere die Méglichkeit von Resonanzen 
zwischen freien und erzwungenen Schwingungen, wobei aber keine definitiven 
Aussagen gemacht werden kénnen und offensichtlich die dynamische Theorie 
herangezogen werden muss. Die entwickelte Theorie 2. Ordnung ist von wesentlicher 
Bedeutung fiir die Theorie enger Doppelsterne und die genaue Erfassung der Bewe- 
gung kiinstlicher Erdsatelliten. Die Theorie sollte beispielsweise in der Lage sein, 
die freie Gestalt der Erdoberflache (als Aquipotentialflache) mit einer Genauigkeit 
von der Gréssenordnung 20 cm wiederzugeben. E, RotH-DESMEULES 
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Grenzschichtforschung - Boundary Layer Research. Symposium der 
Internationalen Union fiir Theoretische und Angewandte Mechanik (IUTAM), 
Freiburg i. Br., 26. bis 29. August 1957. Herausgegeben von H. GORTLER (Springer- 
Verlag, Berlin 1958). 411 S., 206 Abb.; DM 67.50. : } 

Ein Gremium von 90 Strémungs-Wissenschaftlern aus 17 Landern versammelte | 
sich zur ersten IUTAM-Arbeitstagung 1957 in Freiburg i. Br. Wahrend 8 Sitzungen 
gaben 32 Vortragende und 21 Diskussionsredner ihre Beitrage zu den Fortschritten . 
der Grenzschichtforschung bekannt. _ 

Wahrend der Tagung wurde das Interesse auf die physikalischen Probleme der 
Grenzschichtforschung konzentriert, obschon dieser Zweig der Str6mungsmechanik — 
fiir die Technik von massgebender Bedeutung ist. Die inhaltsreichen Manuskripte 
der Referenten und die Diskussionsbeitrage sind in diesem Buch zusammengestellt 
und in der gewohnten, mustergiiltigen Art beim Springer-Verlag gedruckt worden. 
Im Rahmen einer kurzen Buchbesprechung ist es nicht méglich, auf Details des von 
den zahlreichen Autoren behandelten Stoffes einzugehen; kurze Hinweise auf die 
Themen (und Angabe der Verfasser) mégen zur Einfiihrung dieses bedeutenden 
Sammelwerkes gentigen: ; 

Beitrage zur turbulenten Reibungsgrenzschicht (TOWNSEND, BJORGUM, WIEG- 
hardt), spezielle stationare Grenzschichten (ACKERET, LUDWIG, SETH, STEWARTSON, 
GLAUERT, THWAITES, GORTLER), Stabilitatstheorien und Umschlagsbedingungen 
(SCHUBAUER, GORTLER und WiITTING, ZaaT, LAUFER, PFENNINGER, DRYDEN, 
WILLE, ScHoLz, Tatsumi, Lin, Howarp, EICHELBRENNER und MICHEL, BERGH, 
van Driest), Temperatur- und Reibungsgrenzschichten bei Uberschallstromungen 
und in Gasen kleiner Dichte (OSTRACH, SCHUH, PoPPER und REINER, SCHAAF und 
TaLsBot, Lunc und LuBoNSKI, KAPPELER und ZADDACH, SCHERBERG), Wechsel- 
wirkung zwischen Grenzschichten und Verdichtungsstéssen (GADD, FERRARI, 
Mires, MassiGNon, IBERALL), Grenzschichten bei instationaren Vorgangen und 
bei periodischer Absaugung (MoorE, PERSEN, WUEST, PFENNINGER, WORTMANN, | 
Tani), Grenzschichttheorie bei entstehender Ablésung und fiir Zweikomponenten- 
Flissigkeit (NIKOLSKIJ, RACHMATULIN), dreidimensionale Grenzschichten, Ablése- 
bedingungen, laminar-turbulenter Ubergang im Freistrahl und anderes (TrmMan, 
OswatitscH, Watz, TANI, WEHRMANN und WILLE, NICKEL, FROSSLING, DE KRa- 
SINSKI). 

Physiker und Ingenieure, die sich mit den Phanomenen der Reibungsgrenzschich- 
ten zu befassen haben, finden in dem Buche nicht nur neueste Forschungsergebnisse, 
sondern auch die dazu gehérenden Literaturhinweise. Die Beitrage sind in deut- 
scher oder englischer Sprache geschrieben. H. SPRENGER 


Physical Climatology, 2. Aufl. Von Hetmut LanpsBerG, U. S. Weather 
Bureau, Washington (Gray Printing Co., Inc., DuBois, Pennsylvania, 1958). 446 S., 
109 Fig.; $6.00. 

Die erste Auflage dieser grundlegenden Einfiihrung in die Klimatologie erschien 
1941. Nachdem das Werk in kurzer Zeit vier Neudrucke erlebt hatte, wurde es in 
der jetzt vorliegenden zweiten Auflage einer ebenso notwendigen wie griindlichen 
Revision unterzogen. Das Buch, das keine Kenntnisse der héheren Mathematik 
voraussetzt, gliedert sich in vier Kapitel. Im ersten Kapitel, das sich mit der 
Sammlung und dem Gebrauch von klimatologischen Daten befasst, werden neben 
den Untersuchungsmethoden und der Einrichtung der Instrumente, vor allem auch 
die statistischen Methoden behandelt. Das zweite Kapitel geht auf die klimatolo- 
gischen Elemente selber ein. Es ist erfreulich, dass die immer wichtiger werdende 
Warmebilanz der Atmosphare hier eine eigene Darstellung erfahrt. Unter der Ver- 
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bindung der klimatologischen Elemente, die den Inhalt des dritten Kapitels aus- 
macht, versteht H. Lanpspere ihre Vereinigung zu Bildern des Klimas der Welt, 
grosser Regionen der Erde oder auch ganz kleiner Bezirke, also des Klimas der 
‘bodennahen Luftschichten im Sinne von R. GricEr. Das letzte Kapitel schliesslich 
handelt von der angewandten Klimatologie. 

Das sehr anschaulich und instruktiv geschriebene Buch ist nicht nur eine 
ausgezeichnete Einfiihrung in die Klimatologie fiir Studierende der Meteorologie; 
auch der Fachmann wird es mit Vorteil lesen. J.C. THams 


Theorie und Anwendung der direkten Methode von Ljapunov. Von 
WoLrcancG Haun. Ergebnisse der Mathematik und ihrer Grenzgebiete. Neue Folge, 
Heft 22 (Springer-Verlag, Berlin-Géttingen-Heidelberg 1959). 142 S., DM 28.-. 

Die Frage nach der Stabilitat spielt bei Bewegungen und vor allem bei Schwin- 
gungsproblemen, insbesondere in der Regelungstechnik, praktisch eine hervor- 
ragende Rolle. Zu ihrer Untersuchung erweist sich die zweite oder direkte Methode 
von Lyapunov (die die Kenntnis der Lésungen der Differentialgleichungen nicht 
voraussetzt) als sehr geeignet, und sie lasst sich auch zum Aufbau einer allgemeinen 
Stabilitatstheorie heranziehen. Die Theorie der direkten Methode ist in den letzten 
Jahren vor allem durch russische Mathematiker stark geférdert worden, und da 
heute ein gewisser Abschluss erreicht wurde, erscheint eine Ubersicht dringend 
notwendig, besonders auch deshalb, weil manche Originalarbeiten nur schwer zu- 
ganglich sind. Diese Ubersicht vermittelt in ausgezeichneter Weise die Darstellung 
von W. Haun. Die beiden ersten Kapitel entwickeln die elementare Theorie (hin- 
reichende Bedingungen), das dritte bringt Anwendungen auf technische Probleme, 
wahrend die restlichen dem Ausbau der Theorie (Umkehrung der Hauptsatze, 
Empfindlichkeit gegen Sto6rungen, kritische Falle) und dem Aufbau einer allgemei- 
nen Stabilitaétstheorie gewidmet sind. Im Vordergrund stehen naturgemass die 
Stabilitatskriterien fiir gewOhnliche Differentialgleichungen; nur das letzte Ka- 
pitel bringt die Erweiterung auf partielle Differentialgleichungen und Differenzen- 
gleichungen. Die Darstellung schliesst mit einem Literaturverzeichnis, das bis zum 
Jahre 1957 méglichste Vollstandigkeit anstrebt. Das Buch fillt im deutschen 
Sprachgebiet eine empfindliche Liicke aus, und sein Erscheinen ist daher sehr zu 
begriissen. E. RotH-DESMEULES 


Plasma Physics and the Problem of Controlled Thermonuclear Re- 
actions, Vol. III. Herausgegeben von M. A. Leontovicu. Ubersetzung von J. B. 
SyKeEs, Harwell. (Pergamon Press, London 1959). 422 S., 114 Fig. £8.-. 

Der vorliegende Band ist Teil einer vierbandigen Reihe und enthalt eine Serie 
von kiirzlich zur Publikation freigegebenen Arbeiten des Institutes fiir Kernphysik 
der USSR. Die meisten Artikel stammen aus den Jahren 1955-1957 und vermitteln 
ein eindrucksvolles Bild von den Fortschritten Russlands auf dem Gebiet der 
Fusionsforschung bis zum Zeitpunkt der Genfer Konferenz 1958. Der dritte Band 
enthalt 26, vorwiegend theoretische Arbeiten, an denen etwa 20 Autoren beteiligt 
sind. Die behandelten Probleme sind in erster Linie die Erzeugung hoher Tempe- 
raturen, die magnetische Fiihrung geladener Teilchen, Rekombinationseffekte in 
ionisierten Gasen, Strahlungsverluste usw. 

G. I. BupKER diskutiert Methoden zur direkten Umwandlung thermonuklearer 
Energie in elektrische Arbeit. Einen besonders interessanten Beitrag liefert B. A. 
TRUBNIKOW, in dem er zeigt, dass die magnetische Bremsstrahlung einer Plasma- 
schicht die Minimalgrésse eines selbsterhaltenden Fusionsreaktors bestimmt und 
dass diese kritischen Dimensionen selbst im Falle eines D-T—Reaktors sehr betracht- 
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lich sind. Einer der wenigen experimentellen Artikel ist der von N. V. F ILIPPOY, i 
dem er eine Methode zur Messung von Drucken in gepulsten Gasentladungen mi : 
Hilfe von Piezoelementen beschreibt. L. I. RupaKov und R«Z. SADGEEv schlagen 
eine Methode der Plasmaerhitzung durch Hochfrequenz vor, die von der Zyklotron- 
resonanz Gebrauch macht. In einer Arbeit von V. I. Kocan wird die Rekombina- 
tionsstrahlung aus einem Wasserstoffplasma berechnet und mit der Intensitat der 
Bremsstrahlung verglichen. S. I. BRapinski1 und B. B. KapomtseEv zeigen, dass es 
mit Hilfe einer «Guardring»-Methode méglich sein sollte, die Form einer Plasma- 
oberflache zugunsten einer méglichst stabilen Konfiguration zu beeinflussen. 
A. A. VEDENov und R. Z. SADGEEV weisen die Méglichkeit auf, dass fiir eine aniso- 
trope Geschwindigkeitsverteilung Plasmainstabilitaten auftreten kénnen, die bei 
der hydromagnetischen Behandlung in alteren Arbeiten eventuell ibersehen wur- 
den. Im Ganzen beurteilt, diirfte dieses Buch ein neuer Beweis dafiir sein, dass die 
vielseitigen Bemiihungen des Pergamon Institutes, dem westlichen Leser die Lite- 


ratur der USSR zuganglich zu machen, ausserordentlich fruchtbar sind. 
F. HEINRICH 


Kommunikation und Kybernetik in Einzeldarstellungen. Band 1: 
Grundlagen und Anwendungen der Informationstheorie. Von Prof. Dr. W. MEYER- 
EppLer. (Springer-Verlag, Berlin 1959). 446 S., 178 Abb., 1 Tafel; DM 98.-. 

Das Buch ist als erster Band einer Reihe iber Kommunikation und Kyberne- 
tik in Einzeldarstellungen erschienen, deren Herausgeber ebenfalls der Autor war. 
Es ist nur zu hoffen, dass durch den leider zu friih erfolgten Hinschied von 
W. MEYER-EPPLER (Juli 1960) die Herausgabe der Reihe keinen Unterbruch er- 
fahrt und die folgenden Bande eine ahnliche Hohe der Qualitat aufweisen werden. 
Die Darstellung der Informationstheorie hatte in keine bessern Hande gelegt wer- 
den kénnen. Es ist erstaunlich, wie originell und gut fasslich die bereits bestehende 
zahlreiche Literatur mithinein verarbeitet wurde. Der Quellennachweis wurde 
aufs sorgfaltigste geftihrt; und doch erscheint das Buch wie aus einem Guss. 

Inhaltlich ist es in die folgenden elf Kapitel gegliedert, wobei die Klammern 
die Seitenzahlen angeben: Die Kommunikationskette (5) — Strukturtheorie der 
Signale (34) — Eigenschaften lineaver Ubertragungssysteme (14) — Symbolstatistik 
(80) — Gestorte Systeme (20) — Sicherung gegen Ubertragungsfehler (19) — Die Sinnes- 
organe als Informationsempfdanger (78) — Signal und Zeichen (32) — Akustische und 
optische Valenzklassen als Zeichentriiger (47) — Formstrukturen und Konstruktionen 
(25) — Die gestérte sprachliche Kommunikation (48) — Anhang mit einer Tabelle der 
Funktion -pldp und getvenntem Namen- und Sachverzeichnis (15). 

Besonders ausfiihrlich — und mit Recht — sind die Kapitel iiber Symbolstatistik 
und die Sinnesorgane als Informationsempfanger gehalten. Im erstern sind die 
Kernsatze der Informationstheorie am Beispiel der Sprache entwickelt, und im 
andern finden die Eigenschaften der Sinnesorgane beziiglich ihrer Zellenstruktur, 
ihrer Reizschwellen, ihrer Unterscheidbarkeitsgrenzen und der Darstellung der 
Reizempfindungen (Farbdreieck, Phoneme, Vokaldreieck) ihren Platz. Erst 
durch die eingehende Kenntnis beider lasst sich die Kommunikationskette infor- 
mationstheoretisch untersuchen, was in den dem oben genannten Kapitel folgen- 
den durchgefiihrt wird. 

Das vorliegende Buch schliesst eine Liicke im Schrifttum, es befriedigt nicht 
nur den Theoretiker, sondern auch den Fernmeldeingenieur, und selbst der Lin- 
guist wird es mit grossem Gewinn zu seinen Arbeiten verwenden kénnen. 
W. MEYER-EppLeER hat sich mit diesem Buch selbst das schénste Denkmal ge- 
setzt. H. WEBER 
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1 Geomeirische Grundlagen. Einleitung; Koordinatentransformation; Skalar und 
_Vektor; Tensor; Verallgemeinerung; Der e-Tensor; Hauptachsen eines symme- 
_trischen Tensors zweiter Stufe; Tensorfelder; Integralsitze; Krummlinige Koordi- 
-naten..— JJ. Spannungszustand. Spannungstensor; Gleichgewichtsbedingungen ; 
_Hauptspannungen und Hauptschubspannungen; Mohrsche Darstellung von Span- 
“nungszustanden; Ebener Spannungszustand; Cauchysche Spannungsflachen; Hy- 
‘ drostatik. — III. Bewegungszustand. Drehungs- und Verformungsgeschwindigkeit; 
_Ebener Bewegungszustand; Integrabilitatsbedingungen; Materielle Anderungs- 
_geschwindigkeiten. — IV. Grundgesetze. Erhaltung der Masse; Impulssatz; Drallsatz; 
Energiesatz; Einfache Stoffgleichungen. — V.Ideale Flissigkeiten. Bewegungs- 
“gleichungen und Wirbelsdtze; Stationdre Bewegung; Wirbelfreie Stré6mung. — 
VI. Zihe Flissigheiten. Grundgleichungen und Ahnlichkeitsgesetz; Unzusammen- 
“driickbare Fliissigkeit: Strenge Lésungen; Unzusammendriickbare Fliissigkeit: 

Grenzschichtgleichungen; Unzusammendriickbare Fliissigkeit: Nicht-Newtonsches 
Verhalten. — VII. Zéah-plastische und ideal-plastische Stoffe. Stoffgleichungen; Bewe- 
-gungsgleichung des zah-plastischen Stoffs; Grenzschichtgleichungen des zab- 
plastischen Stoffs; Ebener Flie8zustand. — VIII. Hypo-elastische Stoffe. Klassische 
| Elastizitétstheorie. Spannungsgeschwindigkeit; Hypoelastische Stoffe; Hypoelasti- 
‘sche Stoffe in der Nachbarschaft des spannungslosen Zustands; Grundgleichungen 
‘der Elastostatik; Torsion zylindrischer Stabe; Extremalprinzipe; Elastische 
‘Wellen. — IX. Endliche Verzerrung. Almansischer Verzerrungstensor; Greenscher 
‘Verzerrungstensor; Andere Verzerrungstensoren; Lagrangescher und Kirchhoff- 
-scher Spannungstensor. — X. Elastische und hyperelastische Stoffe. Elastische Stoffe; 
‘Hyperelastische Stoffe; Unzusammendriickbare hyperelastische Stoffe; Eindeutig- 
keit und Stabilitat. 


Das Buch soll dem in steigendem Masse gefiihlten Bediirfnis nach einer Einleitung 
in die Kontinuumsmechanik entgegenkommen, die einerseits als Grundlage ftir das 
-vertiefte Studium von Sondergebieten wie Hydrodynamik, Gasdynamik, Elastizitat 
‘und Plastizitat dienen kann und andererseits dem Leser, der kein solches Studium 
plant, eine Einsicht in typische Schlussweisen und Problemstellungen darbietet. 
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Actes du Symposium de Rome. (20-24 septembre 1960) 
Organisé par le Centre International Provisoire de Calcul 


(1960) 680 pp. Fr. 35.— (DM 35,—) 


Authors/Auteurs: M. Altmann, Z. Alterman, H. A. Antosiewicz, N. Artemiadis, 
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